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1.  Statement  of  the  Problem  Studied. 


The  research  was  on  the  effects  of  adsorbing  and  desorbing  molecules  on  the 
mean  resonant  frequency  and  frequency  fluctuations  of  thickness-shear  resonators, 
in  particular,  10  MHz  and  525  MHz  AT-cut  quartz  resonators.  The  study  examines 
the  effects  of  temperature  and  pressure  on  changes  in  mean  resonant  frequency 
and  frequency  fluctuations  due  contaminant  molecules  on  electrode  surfaces  of 
thickness-shear  resonators.  One  or  more  species  and  multilayers  of  adsorbing 
molecules  are  treated. 

The  effect  on  resonant  frequency  due  to  an  adsorbed  thin  mass  layer,  such 
as  one  or  more  layers  of  molecules,  is  measurable  and  is  employed  in  quartz 
microbalances  for  surface  adsorption  studies.  Physical  adsorption  and  desorption 
of  gas  molecules  from  the  resonator  electrodes  are  dependent  on  the  gas  pressure, 
temperature,  heats  of  adsorption  and  sticking  coefficients.  For  a  given  resonator, 
the  mean  resonant  frequency  will  change  as  a  function  of  the  number  of  adsorbed 
molecules,  which  in  turn  is  a  function  of  pressure  and  temperature.  Since  the 
adsorbed  molecules  are  known  to  have  finite  lifetimes  on  the  electrode  surface, 
the  number  of  adsorbed  molecules  at  a  given  instant  will  fluctuate.  Hence,  short¬ 
term  frequency  instabilities  will  result  if  the  resonator  is  sufficiently  sensitive  to 
fluctuations  in  the  number  of  adsorbed  molecules. 

2.  Theoretical  Derivations  and  Details 
of  Calculations. _ 

The  details  of  theoretical  derivations  and  calculations  are  provided  in  technical 
publications  elsewhere.  The  publications  are  listed  in  Section  4. 

The  interested  reader  is  referred  to  the  appendix  for  copies  of 
these  three  papers: 

1.  Derivations  and  calculations  for  the  effects  of  a  monolayer  or  less  of  a  single 
species  of  contaminant  molecules  is  given  in  publication  1. 

2.  Derivations  and  calculations  for  the  effects  of  several  species  of  contaminant 
molecules  in  a  monolayer  or  less  of  contaminant  molecules  is  given  in 
publication  4. 

3.  Derivations  and  calculations  for  the  effects  of  several  layers  of  a  single  species 
of  contaminant  molecules  is  given  in  publication  5. 
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3.  Summary  of  the  Most  Results. 


The  following  conclusions  are  applicable  to  frequency  fluctuations  and 

changes  in  mean  resonant  frequency  induced  by  adsorbing  and  desorbing  con¬ 
taminant  molecules  in  thickness-shear  resonators. 

1.  Calculations  performed  for  a  10  MHz  thickness-shear  resonator  show  that  the 
mean  resonant  frequency  can  change  as  much  as  10  ppm  with  temperature. 
In  practice,  this  frequency  change  is  usually  masked  by  changes  in  material 
properties  of  quartz,  such  as  elastic  stiffnesses  and  density,  which  are  also  tem¬ 
perature  dependent.  The  mean  resonant  frequency  is  also  affected  in  the  same 
order  of  magnitude  by  changes  in  the  gas  pressure.  Frequency-temperature 
and  frequency-pressure  curves  for  the  10  MHz  resonator  cannot  be  described 
by  cubic  polynomial  equations.  (Cubic  polynomial  equations  are  used  in  the 
quartz  resonator  industry  to  describe  frequency-temperature  behavior.)  The 
gas  pressure  would  be  a  more  suitable  parameter  than  temperature  since  the 
quartz  material  properties  is  not  significantly  affected  by  changes  in  pressure. 

2.  Frequency  fluctuations  due  to  adsorbing  and  desorbing  molecules  are  not  a 
simple  function  of  the  percent  area  contaminated.  A  very  “clean”  resonator 
and  a  “dirty”,  namely,  highly  contaminated,  resonator  can  both  have  low 
frequency  fluctuations. 

3.  The  magnitudes  of  frequency  fluctuations  due  to  surface  contaminations  is 
significant  in  VHF,  UHF  and  thin-film  thickness-shear  resonators.  For  a  10 
MHz  resonator,  the  frequency  fluctuations  are  generally  very  small  (<  -160 
dBc  Hz2/Hz)  if  a  monolayer  or  less  of  contaminant  molecules  is  assumed.  If 
the  contaminant  molecules  are  allowed  to  form  multilayers,  calculations  for 
the  10  MHz  resonator  show  that  the  magnitude  of  mean  square  frequency 
fluctuations  at  1  Hz  is  in  the  measurable  range  of  about  -120  dBc  Hz2/Hz. 

4.  The  spectral  density  of  frequency  fluctuations  at  1  Hz  for  multilayer  contam¬ 
ination  is  remarkably  constant  with  changes  in  temperature  if  the  heats  of 
adsorption  for  the  second  and  subsequent  layers  are  similar  to  heat  of  adsorp¬ 
tion  of  the  first  layer.  For  monolayer  adsorption  and  desorption,  the  spectral 
density  of  frequency  fluctuations  at  1  Hz  changes  strongly  with  temperature. 

5.  Magnitudes  of  frequency  fluctuations  caused  by  surface  contamination  can  be 
altered  by: 

a.  changing  the  pressure  of  contaminant  gas 
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b.  changing  the  temperature  and  treating  the  electrode  surface  to  obtain 
different  rates  of  adsorption  and  desorption. 

c.  changing  the  area  of  the  electrodes.  The  smaller  the  electrode  area,  the 
higher  the  frequency  fluctuations. 

6.  If  all  parameters  except  the  plate  thickness  are  held  constant,  the  spectral 
density  of  frequency  fluctuations  is  inversely  proportional  to  the  fourth  power 
of  plate  thickness.  Since  the  resonator  frequency  is  inversely  proportional 
to  plate  thickness,  the  spectral  density  is  also  proportional  to  the  fourth 
power  of  the  resonator  frequency.  Furthermore,  if  the  electrode  diameter 
decreases  linearly  with  the  operating  resonator  frequency,  the  spectral  density 
of  frequency  fluctuations  is  approximately  proportional  to  the  sixth  power  of 
the  resonator  frequency. 

7.  Graphs  of  spectral  density  of  frequency  fluctuations  as  a  function  of  Fourier 
frequency  exhibit  slopes  ranging  from  l/f°  to  1/f2.  1/f  flicker  noise  cat-  be 
obtain  if  the  electrode  surface  has  a  wide  range  of  values  of  heat  of  adsorption 
for  contaminant  molecules. 

4.  List  of  Publications. 
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Kong  Yong  and  John  R.  Vig,  Proceedings  of  the  1989  Ultrasonics  Sympo¬ 
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Abstract— The  mr  ss  loading  effects  of  adsorbing  and  desorbing  con¬ 
taminant  molecules  on  the  magnitude  and  characteristics  of  frequency 
fluctuations  ir  a  thickness-shear  resonator  are  studied.  The  study  is 
motivated  by  the  observation  that  the  frequency  of  a  thickness-shear 
resonator  is  determined  predominantly  by  such  mechanical  parame¬ 
ters  as  the  thickness  of  the  resonator,  elastic  stiffnesses,  mass  loading 
of  the  electrodes,  and  energy  trapping.  An  equation  was  derived  relat¬ 
ing  the  spectral  density  of  frequency  fluctuations  to:  1)  rates  of  ad¬ 
sorption  and  desorption  of  one  species  of  contaminant  molecules,  2) 
mass  per  unit  area  of  a  monolayer  of  molecules,  3)  frequency  constant, 
4)  thickness  of  resonator,  and  5)  number  of  molecular  sites  on  one  res¬ 
onator  surface.  The  induced  phase  noises  were  found  to  be  significant 
in  very-high-frequency  resonators  and  are  not  simple  functions  of  the 
percent  area  contaminated.  The  spectral  density  of  frequency  fluctua¬ 
tions  was  inversely  proportional  to  the  fourth  power  of  the  thickness  if 
other  parameters  were  held  constant.  Since  the  resonator  frequency  is 
inversely  proportional  to  the  thickness,  the  spectral  density  is,  in  effect, 
proportional  to  the  fourth  power  of  resonator  frequency. 


I.  Introduction 

SHORT-TERM  FREQUENCY  STABILITY  (i.e., 
noise)  is  among  the  least  well-understood  resonator 
phenomena.  Factors  that  have  been  recognized  as  sources 
of  short-term  instabilities  include:  Johnson  noise,  tem¬ 
perature  fluctuations,  random  vibrations,  acoustic  losses, 
noise  originating  from  interfaces  between  the  electrodes 
and  the  quartz  plate  and  between  the  mounting  structure 
and  the  quartz  plate,  and  noise  due  to  the  oscillator  cir¬ 
cuitry.  Vig  [1]  proposed  a  prelimimuy  model  of  frequency 
fluctuations  induced  by  surface  contaminations  and  de¬ 
scribed  qualitatively  the  effects  of  such  contaminations. 

This  paper  explores  and  estimates  the  frequency  fluc¬ 
tuations  caused  by  adsorption  and  desorption  of  one  spe¬ 
cies  of  molecules  on  the  surface  of  a  thickness-shear 
quartz  resonator.  The  characteristics  of  such  frequency 
fluctuations  are  studied.  Our  study  is  motivated  by  the 
observation  that  the  resonant  frequency  of  a  thickness 
shear  resonator  is  determined  predominantly  by  such  me¬ 
chanical  parameters  as  the  thickness  of  the  resonator, 
elastic  stiffnesses,  mass  loading  of  the  electrodes,  and  en- 
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ergy  trapping.  Hence,  the  effects  on  the  mass  loading  of 
the  stochastic  characteristics  of  mass  loading  due  to  ad¬ 
sorption  and  desorption  of  contaminant  molecules  are 
considered.  These  effects  are  confined  to  the  electrode 
area,  since  the  vibrational  energy  is  predominantly 
trapped  within  the  electrodes. 

II.  Adsorption/Desorption  of  Contaminant 
Molecules 

A  dynamt  equilibrium  exists  between  the  resonator 
surface  and  the  gas  above  it.  This  equilibrium  is  main¬ 
tained  by  the  rate  of  adsorption  and  desorption  of  contam¬ 
inant  molecules  on  the  resonator  surface.  The  area  of  this 
surface  is  taken  as  the  electrode  area,  since  most  of  the 
resonator  energy  is  confined  within  the  electrodes  (energy 
trapping). 

The  rate  of  arrival  of  molecules  at  the  surface  can  be 
evaluated  readily  from  simple  kinetic  theory  of  gases.  A 
simple  form  of  the  relationship  between  the  rate  of  ar¬ 
rival,  pressure,  temperature  and  molecular  weight  of  the 
molecule  [2]  is 

r  =  3.51  x  1022  P /{TM)'^1  molecules/cm2/s  (1) 

where  P  is  the  pressure  in  torrs,  T  the  temperature  is  Kel¬ 
vins,  and  M  the  molecular  weight.  In  the  enclosed  space 
of  a  crystal  resonator,  the  pressure  will  vary  slightly  with 
the  fluctuations  in  the  number  of  adsorbed  molecules.  For 
the  present  linear  model,  a  constant  pressure  is  assumed. 
If  we  consider  the  following  example: 

P  =  1  torr 

T  =  300°K 

M  =  28  (nitrogen), 

the  rate  of  arrival  r0  of  nitrogen  molecules  at  a  contami¬ 
nant  site  with  a  spacing  of  0.5  nm  is  r0  -  9.6  >  105 
molecules/site/s.  For  a  given  species  of  molecules  at  a 
fixed  temperature,  the  rate  would  be  linearly  dependent 
on  the  pressure.  The  reciprocal  of  this  rate  is  the  mean 
interarrival  time  r0. 

The  contaminant  moiecures  reside  on  the  surface  for  a 
finite  time.  In  a  simple  physical  adsorption  process,  the 
mean  residence  time  13}  is 

ri  =  ct  exp  ( Ej/(RT ))  (2) 

where  R  is  the  gas  constant  =  1.987  x  10'3  keal/mol  K, 
T  is  temperature  in  Kelvins,  and  Ed  is  the  desorption  en- 
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ergy  in  kcal/mol.  The  exact  value  of  the  parameter  a  is 
conceptually  not  essential  to  our  derivations.  For  subse¬ 
quent  calculation  purposes,  we  assign  to  a  a  value  of  t  x 
10" 13  seconds.  If  we  use  a  typical  range  of  values  for  Ed 
of  10  to  20  kcal/mol,  the  mean  residence  time  is  rt  =  2 
x  10"6  to  50  s.  The  mean  desorption  rate  of  molecule 
from  the  surface  r,  is  the  reciprocal  of  this  mean  residence 
time. 

A  monolayer,  or  less,  of  adsorbed  contaminant  mole¬ 
cules  is  assumed.  When  the  system  of  contaminant  gas 
and  resonator  surface  settles  down  to  a  steady-state  dy¬ 
namic  equilibrium,  the  states  of  the  contaminant  sites, 
which  are  either  unconf"*minated  (state  0)  or  contaminated 
(state  1),  can  be  described  by  steady-state  probabilities. 
A  transition  diagram  in  Fig.  1  shows  the  interaction  of  the 
two  states  in  a  site.  The  rate  in  which  the  site  enters  a 
particular  state  must  be  equal  to  the  rate  at  which  it  leaves. 
Hence,  balance  equations  for  the  steady-state  probabili¬ 
ties  can  be  written  as  follows: 


State 

Rate  Entering 

Rate  Leaving 

0 

P\r\ 

=  Por0 

1 

Poro 

=  Pin 

where  p0  and  p,  are  the  steady-state  probabilities  of  state 
0  and  state  I,  respectively.  In  writing  these  balance  equa¬ 
tions,  the  assumption  was  made  that  the  molecules  arriv¬ 
ing  at  a  site  would  be  adsorbed  immediately  if  the  site 
were  uncontaminated  but  would  otherwise  bounce  off  the 
surface  and  return  to  the  gas,  that  is,  we  assume  a  sticking 
probability  of  one  if  the  site  is  uncontaminated  and  a 
sticking  probability  of  zero  if  the  site  is  contaminated.  In 
reality,  both  probabilities  generally  have  a  value  between 
zero  and  one.  The  exact  value  of  the  sticking  probability 
will,  depending  on  the  desorption  rate,  lower  or  increase 
the  magnitude  of  frequency  fluctuations,  but  will  other¬ 
wise  not  affect  the  logic  of  our  derivations. 

Since  the  site  must  at  any  given  time  be  in  either  state 
0  or  1 ,  p0  and  pi  must  sum  to  one: 

Po+Pt  =  l.  (4) 


Ill.  Frequency  Change  Due  to  a  Very  Small 
Mass  Loading 

A.  Frequency  Influence  Curve  of  an  SC-Cut  Strip 

The  effect  of  a  very  small  mass  on  the  surface  of  a 
thickness-shear  resonator  is  investigated;  namely,  the  ef¬ 
fect  on  the  resonant  thickness-shear  frequency.  For  this 
purpose  a  one-dimensional  finite-element  model  is  em¬ 
ployed  |4j.  Since  the  ratio  of  the  mass  loading  per  unit 
area  to  the  resonator  mass  per  unit  area  is  expected  to  be 
very  small  (<1  x  10-4),  the  magnitude  of  frequency 
change  is  linearly  proportional  to  the  mass  loading.  Con¬ 
sequently,  the  finite  clement  calculations  for  frequency 
change  at  a  mass  loading  ratio  of,  say,  1  x  10‘4  can  be 
scaled  linearly  to  yield  frequency  change  due  to  mass 
loading  of  contaminant  molecules. 

Fig.  2  shows  a  unit  width  SC-cut  strip  subjected  to  a 
small  mass  loading  (length  99  pm,  thickness  10  nm,  and 
density  equal  to  that  of  quartz).  The  position  of  the  mass 
is  given  by  the  coordinate  along  the  length  of  the  strip. 
The  strip’s  resonant  frequency  is  3.6  MHz.  If  the  fre¬ 
quency  change  can  be  plotted  against  the  position  of  the 
mass  on  the  strip,  a  curve  known  as  the  frequency  influ¬ 
ence  curve  is  obtained.  Fig.  3  gives  the  frequency  influ¬ 
ence  curve  for  a  mass  loading  (thickness  0.5  nm  and 
length  0.5  nm)  obtained  from  scaling  the  finite  element 
calculations  for  the  strip  of  Fig.  2.  We  observe  that  the 
mass  loading  has  predominant  frequency  influence  over 
the  middle  third  portion  of  the  strip.  The  area  under  this 
frequency  influence  curve  represents  the  frequency  change 
due  to  a  monolayer  of  molecules  of  0.5-nm  thickness. 

B.  Frequency  Cha.'.ge  per  Site  over  a  2-D  Surface 

For  simplicity,  we  assume  a  flat  frequency  influence 
surface  over  the  entire  electrode  area,  although  the  actual 
frequency  influence  of  contaminant  molecules  is  predom¬ 
inantly  over  the  middle  third  area  of  electrodes.  Subse¬ 
quently,  it  will  be  shown  that  the  noise  levels  are  in¬ 
versely  proportional  to  the  surface  area  if  other  parameters 
are  kept  constant.  Therefore,  the  assumption  of  a  flat  fre¬ 
quency  influence  surface  will  lower  the  noise  levels.  Based 
on  the  one-dimensional  frequency  influence  curve  of  Sec¬ 
tion  III-A,  this  assumption  resulted  in  noise  level  esti¬ 
mates  that  were  at  least  3  dBc  lower  than  when  the  actual 
frequency  influence  surface  was  used. 

For  a  flat  frequency  influence  surface,  the  change  in 
frequency  due  to  a  monolayer  of  contaminant  molecules 
is 


Therefore,  from  (3)  and  (4)  we  obtain 


Po  - 


— Dll. . 
fo  +  rx 


and 


Pi  = 


r0  +  r,  ‘ 


(5) 


The  probability  p\  can  be  interpreted  as  the  ratio  of  con¬ 
taminated  area  to  total  surface  area,  that  is,  the  fraction 
of  the  area  contaminated.  A  highly  contaminated,  or 
"dirty,”  surface  would  have  a  value  of px  close  to  unity. 


a  f  = 

— J 


f 

JO 


m' 

m 


(6) 


where/,,,  m',  and  m  are,  respectively,  the  resonator  fre¬ 
quency.  mass  per  unit  area  of  a  monolayer  of  molecules, 
and  mass  per  unit  area  of  resonator.  The  frequency  change 
per  site  is  constant 


Si  = 


A/ 

N’ 


i  =  1,2, 


,  N 


(7) 
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VERY  m  mss  LOADING 


Fig.  2.  SC  cut  strip  {resonant  frequency  =  3.6  MHz)  subjected  to  very 
small  mass  loading. 


r reeveocy  Influence  Curve 


o<, 


«0  I.)  it  JO  40  JO  4  0  70  40  40  10  0 

<*o 


Fig.  3.  Frequency  influence  curve  of  SC-cut  strip  (resonant  frequency  = 
3.6  MHz)  of  Fig.  2  due  to  one  contaminated  molecular  site  with  dimen¬ 
sions  0.5  x  0.5  nm,  and  density  same  as  quartz. 


where  N  is  the  number  of  sites  in  the  active  area  of  the 
resonator;  that  is,  the  area  of  electrode.  The  frequency 
change  per  site  is  calculated  for  two  resonators. 


mean  or  expected  value  of  b(  (f)  is  [5]  simply 

E[bi(t)]=pl,  (9) 

where  p\  is  the  probability  of  the  site  being  contaminated 
which,  from  (5),  is  the  same  as  the  fraction  of  area  con¬ 
taminated.  The  variance  of  (t)  is  [5] 

var[M/)]  =  E[b,(t)2]  -  £[M0f 

=  Pi(l  —  Pi)-  (10) 

If  pi  was  substituted  using  (5),  (10)  could  be  expressed  in 
terms  of  the  adsorption  and  desorption  rates  as 

var[M0j  =  2.  (11) 

(>o  +  n) 

A.  Frequency  Fluctuations 

The  resonator  frequency  fluctuations  induced  by  ad¬ 
sorption  and  desorption  of  contaminant  molecules  are  ob¬ 
tained  by  summing  over  2 N  sites  the  stochastic  processes 
fO) 

2N 

/(r)  =2/(0  (12) 

i  ■=  i 

where  the  factor  2  is  included  because  the  thickness  shear 
resonator  has  two  surfaces.  The  term  on  the  right-hand 
side  is  substituted  by  (8)  to  yield 

2N 

/(<)-S|A(i).  (13) 

i  “  I 


Example  1:  10~MHz  Thickness-Shear  Resonator 
Electrode  area  =  10  mm2 
Thickness  =0.165  mm 

0.5-nm  monolayer  of  molecules  (density  =  quartz) 
Number  of  sites  =  4  x  1013  (0.5  nm  spacing) 

A /  =  30  Hz  or  3  ppm 
=  7.5  X  10" 13  Hz 

Example  2:  525-MHz  AT-Cut  Resonator 
Electrode  area  =  0.003  mm2 
Thickness  =  3.2  pm 

0.5-nm  monolayer  of  molecules  (density  =  quartz) 
Number  of  sites  =  12  x  109  (0.5  nm  spacing) 

A /  -  82  kHz  or  156  ppm 
=  6.83  x  10'6  Hz 

We  observe  that  the  frequency  change  per  site  can  in¬ 
crease  by  orders  of  magnitude  if  the  thickness  and  elec¬ 
trode  area  of  resonator  are  reduced. 

IV.  Stochastic  Processes 

Each  contaminant  site  yields  a  stochastic  process  that 
is  assumed  stationary,  and  mutually  independent  of  other 
sites.  The  stochastic  process  of  frequency  fluctuations  per 
site  /  can  be  written  as 

fi(')  =  ZMt)  /  =  1,2,  •  •  •  ,  N,  (8) 


Since  bt(t)  is  assumed  stationary,  and  mutually  inde¬ 
pendent,  the  mean  and  variance  of  /(f)  can  be  derived 
easily  by  using  elementary  principles  of  statistics.  The 
mean  of /(f)  is,  employing  (9), 

*1/(0]  =  2  *f*[M0] 

i  “  1 

=  Px  2  h.  (14) 

t «  1 

Equation  (14)  can  be  written  in  a  simple  form,  since  the 
sum  of  frequency  change  per  site  over  N  sites  is  equal  to 
A/ 

E[f(t)\  =  2p,A/  (15) 

Hence,  the  mean  frequency  fluctuation  is  twice  the  prod¬ 
uct  of  the  fraction  of  area  contaminated  and  frequency 
change  due  to  a  monolayer  of  molecules. 

The  variance  of /(f),  o2,  can  be  derived  similarly, 
using  (10), 

2  N 

va r[/(f)]  =  .2  £2  var[/?i(f)] 

2N 

O2  -  PiO  -  Pi)  2  $?,  (16) 

i=  i 


where  h,  ( f )  is  a  Bernoulli  random  variable  with  a  value  An  expression  in  terms  of  the  rates  of  adsorption  and  de- 
of  one  if  the  site  is  contaminated  and  zero  otherwise.  The  sorption,  frequency  change  of  a  monolayer  of  molecules, 
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and  total  number  of  sites  is  obtained  by  substituting  (11) 
and  (7)  into  the  right-hand  side  of  (16) 

-  7^-2  Wf/N.  (17) 

(n>  +  n) 

Equation  (17)  is  valid  only  if  a  flat  frequency  influence 
surface  is  assumed.  The  variance  is  inversely  proportional 
in  N  if  other  parameters  were  kept  constant. 

B.  Autocorrelation  and  Spectral  Density  Functions 
The  autocorrelation  function  fl(r)  is  defined  [6]  as 

*(t)  =  £[/('  +  r)/(0]  (18) 

where  the  mean  of  / (t)  has  been  removed  and  r  is  the 
time  lag.  Using  (13)  and  assuming  that  the  sites  are  mu¬ 
tually  independent,  (18)  becomes 

2N 

*(t)  =  S  +  t)  MO].  (‘9) 

In  the  steady-state  dynamic  equilibrium  of  one  species 
of  contaminant  molecule  and  one  type  of  site  on  the  res¬ 
onator  surface,  there  is  one  adsorption  rate  r0  and  one  de¬ 
sorption  rate  rt.  Hence,  the  stochastic  process  has  only 
one  correlation  time  rc 

rc  =  ('ll  +  n)"1.  (20) 

The  autocorrelation  function  in  (19)  yields  an  exponential 
autocorrelation  function  noise  [6],  [7] 

2N 

R(t)  =  2  $?  Var[6l(f)]eXp(-|T|/rc).  (21) 

From  (16)  and  (21),  we  obtain 

R{t)  =  cr2  exp  (-|t|/tc).  (22) 

The  one-sided  spectral  density  function  is  de¬ 

rived  by  taking  the  Fourier  transform  [6]  (Wicner- 
Khintchine  relation)  of  (22) 

W)  -  7-^77  (23) 

1  +  (2t vfre) 

where /is  the  Fourier  frequency.  On  the  right-hand  side, 
o 2  is  substituted  for  (17)  to  yield  the  expression 

c  8r0r,(A ff/N 

A/  (n>  +  r,)3  +  47r2/2(r0  +  r,)' 

It  is  worthwhile  to  describe  in  words  the  meaning  of 
(24).  The  equation  states  that  the  spectral  density  of  fre¬ 
quency  fluctuations  caused  by  one  species  of  comamina.it 
molecule  and  site  is:  1)  a  function  of  some  rational  frac¬ 
tion  of  the  rates  of  adsorption  and  desorption,  2)  propor¬ 
tional  to  the  square  of  the  frequency  change  induced  by  a 
monolayer  of  molecules,  and  3)  inversely  proportional  to 
the  number  of  sites.  The  spectral  density  function  exhibits 
a  white  frequency  noise  process  when  /  «  (r0  + 
T|)/(2t)  and  random  walk  in  frequency  when  /  »  (r0 
+  r|).  There  is  a  simple  relation  between  the  spectra! 


densities  of  frequency  fluctuations  and  phase  noise  of  an 
oscillator  using  the  resonator  [8],  [9];  namely, 

S*(f)  =  S„(f)/f\ 

8r0r,(A/)2/A/  1 

(ro  +  r,)J  +  4ff2/2(r0  +  r,)  f2 

The  spectral  density  of  phase  fluctuations  is  more  com¬ 
monly  used  in  the  field. 

V.  Results 

The  spectral  densities  of  phase  fluctuations  for  the  two 
resonators  in  Examples  1  and  2  of  Section  III-B  are  stud¬ 
ied  over  a  wide  range  of  r0  and  r,  values.  Table  I  shows 
the  calculations  for  the  phase  noise  of  an  oscillator  using 
the  two  resonators  at  300°  K.  The  second  column  of  the 
table  was  calculated  using  (1)  and  the  pressure  values  of 
column  1 .  A  molecular  weight  of  28  and  site  spacing  of 
0.5  nm  were  assumed.  Column  4  was  computed  by  em¬ 
ploying  the  desorption  energies  in  column  3  and  the  re¬ 
ciprocal  of  (2).  Columns  5  and  6  are  obtained  from  (5) 
and  (20),  respectively.  The  last  two  columns  are  the  phase 
noises  at  1  Hz,  which  are  obtained  by  employing  (25), 
parameters  from  Examples  1  and  2  of  Section  III-B,  and 
rates  from  columns  2  and  4. 

We  observe  from  the  second  to  last  column  of  Table  I 
that  the  10  MHz  thickness  shear  resonator  yields  phase 
noise  levels  which  are  barely  measurable  by  today’s  ex¬ 
isting  noise  measurement  technology.  The  last  column, 
on  the  other  hand,  shows  that  the  525-MHz  resonator  can 
produce  significant  levels  of  phase  noise.  We  note  that 
different  noise  levels  are  obtained  for  the  same  fraction  of 
area  contaminated. 

In  order  to  study  the  effects  of  fractional  area  contam¬ 
inated  on  phase  noise  levels,  we  employ  (5)  to  rewrite 
(25)  in  the  following  form 


S»(  1)  = 


8T|£iQ  ~  P,)V/A/ 

r]  +  [2ir(l  -  pt)f 


0  <  pi  <  1 


where  the  spectral  density  is  evaluated  at  1  Hz.  If  we  fix 
the  parameters  A/ and  N,  we  can  plot  the  spectral  density 
againt  the  fraction  of  area  the  contaminated  for  a  given 
value  of  rb  i.e.,  surface  desorption  energy.  Fig.  4  shows 
such  a  plot  for  different  values  of  desorption  energy.  The 
parameters  for  the  525  MHz  resonator  from  Example  2  of 
Section  III-B  were  used.  We  observe  that  for  p{  less  than 
0. 1 ,  the  phase  noise  decreases  linearly  with  decreasing  p,. 
In  general,  the  noise  level  decreases  rapidly  when  /?,  is 
greater  than  0.9.  Hence  the  resonator  is  noisiest  when  the 
fraction  of  the  area  contaminated  is  greater  than  0.1  but 
less  than  0.9.  The  magnitude  of  the  peak  noise  level  de¬ 
pends  on  the  value  of  the  desorption  energy,  the  highest 
of  which  occurs  in  the  range  of  Ed  value,  from  17  to  18 
kcal/mol  (the  parameter  a  of  (2)  being  kept  constant  at  1 
x  10‘ 13  seconds).  The  value  of  pt  at  the  peak  noise  level 
approaches  unity  when  Ed  is  greater  than  20  kcal/mol. 

The  effects  of  the  two  state  variables,  the  pressure  and 
temperature,  on  phase  noise  levels  can  be  studied.  Fig.  5 


456 


IEEE  TRANSACTIONS  ON  ULTRASONICS.  FERROELECTRICS:  AND  FREQUENCY  CONTROL.  VOL.  36.  NO.  4.  JULY  1989 


TABLE  I 

Phase  Fluctuations  Due  to  Surface  Contaminations*  in  IO-MHz  and  525-MHz  Thickness-Shear 

Resonators 


10-MHz 

Resonator'1 

525-MHz 

Resonator* 

Pressure 

torr 

FO 

1/s 

■  Ed 
keal/mol 

F| 

1/s 

Fraction 
of  Area 
Contamination 

*Y 

s 

10  log  5,(1) 
rad'/Hz 

10  log  5,(1) 
rad'/Hz 

10’ 

10’ 

5.49 

10’ 

0.500 

5  x  lO-10 

-196 

-93 

10° 

106 

9.61 

106 

0.500 

5  x  10-7 

-166 

-63 

10-3 

103 

13.7 

103 

0.500 

5  x  10'4 

-136 

-33 

io-J 

10 

16.5 

10 

0.500 

5  X  10_J 

-117 

-13 

!<T6 

1 

17.8 

1 

0.500 

0.5 

-117 

-13 

10"7 

io-‘ 

19.2 

io-‘ 

0.500 

5 

-126 

-22 

io-7 

10'3 

22.0 

10"3 

0.500 

500 

-147 

-42 

*A  heavier  contaminant  molecule  will  raise  the  phase-noise  levels. 
‘Resonators  from  Examples  1  and  2  of  Section  Ill-B. 


Fig.  4.  Effects  of  fraction  of  area  contaminated  on  spectral  density  of  phase 
fluctuations  at  1  Hz  (525-MHz  resonator). 


Fig.  5.  Phase  noise-temperature  curves  for  52j-mHz  resonator. 


exhibits  the  phase  noise-temperature  curves  of  the  525- 
MHz  resonator  at  a  pressure  of  0. 1  torr.  We  observe  that 
when  the  temperature  is  increased  from  a  reference  tem¬ 
perature  of,  say,  25°C,  the  phase  noise  levels  can  in¬ 
crease  or  decrease  depending  on  the  values  of  Ed.  The 
turnover  temperatures  of  the  curves  increase  with  increas¬ 
ing  values  of  Ed. 


Fig.  6  shows  the  phase  noise-pressure  curves  of  the 
525  -MHz  resonator  at  a  temperature  of  300°  K.  The  over¬ 
all  shapes  of  the  curves  look  similar  to  the  curves  in  Fig. 
5  but  the  turnover  pressures  decrease  with  increasing  val¬ 
ues  of  Ed.  We  observe  that  a  decrease  in  pressure  does  not 
necessarily  imply  a  decrease  in  noise  levels. 

Fig.  7  shows  the  spectral  density  of  phase  fluctuations 
for  the  525-MHz  resonator  using  (25)  and  r0  and  r,  equal 
to  100  and  0.1  s_1.  For  Fourier  frequencies  of  less  than 
10  Hz,  the  white  frequency  fluctuations  predominate.  The 
random  walk  frequency  noise  is  dominant  at  Fourier  fre¬ 
quencies  greater  than  50  Hz. 

We  can  make  another  important  observation  based  on 
the  parameters  in  (24)  or  (25).  For  this  purpose,  we  re¬ 
place  the  term  A/ in  (24)  using  (6) 


M/)  = 


8f0ri(/om'/m)2//V 
(fo  +  n)3  +  4r2/2(r0  +  r,) 


(27) 


Upon  a  first  glance,  it  would  appear  that  the  spectral  den¬ 
sity  is  proportional  to  the  square  of  the  resonator  fre¬ 
quency  f0,  if  other  parameters  were  held  constant.  But  the 
frequency  of  a  thickness-shear  resonator  is  inversely  pro¬ 
portional  to  the  thickness  of  the  plate 

fo  =  k/h  (28) 


where  h  is  the  resonator  thickness  and  k  is  the  frequency- 
constant  determined  by  the  elastic  constants  and  overtone 
of  thickness  shear  vibration.  The  mass  per  unit  area  is  also 
related  to  the  thickness 


m  =  ph 


(29) 


where  p  is  the  resonator  density.  Substituting  (28)  and 
(29)  into  (27),  we  obtain 


M/)  = 


8r0r,(W)7(M»4) 

(r0  +  r,)3  +  4  x2/2(r0  +  r,) 


(30) 


Hence,  the  spectral  density  is  inversely  proportional  to 
the  fourth  power  of  h.  Since  the  resonator  frequency  is 
inversely  proportional  to  h,  (30)  shows,  in  effect,  that  the 
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Fig.  6.  Phase  noise-pressure  curves  for  525-MHz  resonator. 


Fig.  7.  Spectral  density  of  phase  fluctuations  for  525-MHz  resonator. 

spectral  density  is  proportional  to  the  fourth  power  of  res¬ 
onator  frequency.  Experimental  results  relating  the  flicker 
frequency  fluctuations  to  the  fourth  power  of  quartz  res¬ 
onator  frequency  were  reported  by  Parker  [9]  and  Kroupa 
[10]. 

In  an  actual  thickness  shear  resonator,  the  electrode  area 
usually  decreases  with  increasing  resonator  frequency, 
which  leads  to  decreasing  values  of  //and  increasing  mag¬ 
nitudes  of  phase  noise.  If  the  electrode  diameter  decreases 
with  frequency  by  a  linear  factor,  the  number  of  sites  N 
will  decrease  by  the  square  of  this  linear  factor,  and  the 
spectral  density  of  frequency  fluctuations,  according  to 
(30)  will  be  proportional  to  the  sixth  power  of  the  reso¬ 
nator  frequency.  This  tends  to  agree  with  the  phenome¬ 
nological  law  proposed  by  Gagnepain,  Uebersfeld,  Gou¬ 
jon,  and  Handel  [11];  namely,  that  the  spectral  density  of 
frequency  fluctuations  is  proportional  to  the  reciprocal  of 
the  fourth  power  of  resonator  unloaded  Q-facior 

If  the  product  Quf0  is  constant,  as  is  usually  the  case,  then 
their  spectral  density  is  proportional  to  the  sixth  power  of 
resonator  frequency 


In  a  recent  review  paper,  Weissman  [12]  reported  that 
a  superposition  of  a  number  of  Lorentzian  spectra  similar 
to  (30)  will  yield  flicker  frequency  fluctuations.  When 
more  than  one  species  of  contaminant  molecules  and  sites 
are  present  in  the  contaminant  gas-resonator  surface  sys¬ 
tem,  the  stochastic  process  has  more  than  one  correlation 
time.  Hence,  depending  on  the  distribution  of  the  corre¬ 
lation  times,  the  basic  model  in  (30)  can  produce  flicker 
frequency  fluctuations  over  some  range  of  Fourier  fre¬ 
quencies.  Most  observed  noise  at  1  Hz  is  flicker  fre¬ 
quency  rather  than  white  or  random  walk. 

VI.  Conclusion 

The  following  conclusions  are  applicable  to  phase  noise 
in  thickness  shear  resonators  induced  by  one  species  of 
contaminant  molecule  and  site. 

1)  The  phase  noise  is  not  a  simple  function  of  the  per¬ 
cent  area  contaminated,  as  shown  in  (26).  A  very 
“clean”  resonator,  and  a  “dirty,”  that  is,  highly 
contaminated,  resonator  can  both  have  low  noise 
levels. 

2)  The  phase  noise  due  to  surface  contaminations  can 
be  significant  in  VHF,  UHF,  and  membrane  reso¬ 
nators.  For  the  10-MHz  resonator  only  the  worst 
case  calculations  begin  to  approach  the  observed 
noise  levels.  In  contrast,  for  the  525-MHz  resona¬ 
tor,  the  worst  case  calculations  are  too  high  for  ob¬ 
served  noise  levels. 

3)  Surface-contamination-caused  noise  levels  can  be 
altered  by: 

a)  changing  the  pressure  of  the  contaminant  gas 
above  the  surface, 

b)  changing  the  temperature  and  treating  the  surface 
to  obtain  different  rates  of  adsorption  and  de¬ 
sorption, 

c)  changing  the  area  of  the  electrode.  The  smaller 
the  area,  the  higher  the  noise  levels. 

4)  If  all  parameters  except  the  thickness  h  are  held  con¬ 
stant  in  (30),  the  spectral  density  of  frequency  fluc¬ 
tuations  is  inversely  proportional  to  the  fourth  power 
of  h.  Since  the  resonator  frequency  is  inversely  pro¬ 
portional  to  h,  the  spectral  density  is  also  propor¬ 
tional  to  the  fourth  power  of  the  resonator  fre¬ 
quency.  Furthermore,  if  the  electrode  diameter 
decreases  linearly  with  the  operating  resonator  fre¬ 
quency,  the  spectral  density  of  frequency  fluctua¬ 
tions  is  approximately  proportional  to  the  sixth 
power  of  the  resonator  frequency. 

One  possible  experiment  to  verify  the  findings  in  this  pa 
per  may  involve  measuring  the  phase  noise  of  a  ultra-high- 
frequency  resonator  while  it  is  in  a  controlled  atmosphere 
in  a  vacuum  chamber,  and  then  decreasing  the  pressure 
gradually.  Extra  care  must  be  taken  to  keep  constant  all 
other  parameters  in  (30)  and  to  reduce  other  noise  sources, 
such  as  random  vibrations,  temperature  fluctuations,  noise 
due  to  oscillator  circuitry,  etc.,  so  as  not  to  interfere  with 
the  measurements.  The  gas  above  the  resonator  must  con- 
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sist  of  essentially  one  species  of  contaminant  molecules. 
A  concomitant  experiment  would  entail  keeping  the  pres¬ 
sure  constant  and  varying  the  temperature  instead.  The 
effects  of  decreasing  thickness  on  phase  noise  levels  in  a 
given  high-frequency  resonator  should  also  be  studied. 
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Abstract 


Resonator  frequency  fluctuations  due  to  adsorption  and  desorption  of 
molecules  on  plate  electrodes  are  studied  using  the  principle,  of  mass-loading 
effects  of  adsorbed  molecules.  The  study  is  based  on  a  525  MHz,  AT-cut  quartz 
resonator  enclosed  in  a  small  crystal  holder.  The  maximum  root  mean  square 
of  pressure  fluctuations  at  300  K  in  a  crystal  holder  with  a  height  and  diameter 
of  respectively  1  mm  and  2mm  is  estimated  to  be  in  the  order  10-6  torr,  and 
hence  would  be  a  factor  in  frequency  fluctuations  if  the  crystal  holder  pressure  is 
of  the  same  order  of  magnitude.  Equations  relating  the  surface  adsorption  rates 
of  crystal  holder  to  pressure  are  derived,  and  found  to  be  quadratic  polynomial 
functions  of  the  adsorption  rates.  Calculations  based  on  these  equations  show 
that  a  contaminant  gas  with  a  higher  desorption  energy  creates  larger  changes  in 
pressure  when  the  temperature  is  varied.  The  function  describing  ti  e  frequency 
fluctuations  due  to  any  one  contaminant  site  is  a  continuous-time  Markoff  chain. 
Kolmogoroff  equations  and  an  autocorrelation  function  for  the  Markoff  chain  are 
derived.  The  autocorrelation,  and  spectral  density  function  of  resonator  frequency 
fluctuations  are  derived.  The  spectral  density  of  frequency  fluctuations  at  1  Hz  is 
studied  as  a  function  of  pressure,  temperature,  and  desorption  energy  of  molecules. 
The  noise  levels  for  a  contaminant  gas  with  one  type  of  molecules  are  found  to 
be  lower  for  lower  desorption  energies,  and  higher  at  lower  pressures.  Graphs 
of  the  power  spectral  density  functions  yield  power-law  noise  processes  which 
range  from  l/f°  to  1/f2.  The  noise  magnitude  is  sensitive  to  the  composition  of 
the  contaminant  gas. 
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I.  Introduction. 


The  effect  of  a  very  thin  mass  layer,  such  as  a  monolayer  of  molecules,  on  the 
frequency  of  a  quartz  resonator  is  measurable.  This  principle  is  employed  in  quartz 
microbalances  for  surface  adsorption  studies[l].  Since  the  molecules  adsorbed  on 
the  surface  am  known  to  have  finite  lifetimes,  the  number  of  adsorbed  molecules 
at  a  given  instant  will  fluctuate.  Hence,  short-term  frequency  instabilities  will 
result  if  the  resonator  is  sufficiently  sensitive  to  a  monolayer  of  molecules. 
This  was  shown  to  be  theoretically  possible  in  ultra-high  frequency  thickness- 
shear  resonators[2].  The  effect  may  be  greater  in  thin-film  resonators  where  the 
thickness  of  the  vibrating  element,  and  the  size  of  electrode  patch  are  critical. 
This  paper  examines  the  magnitude  and  spectral  characteristics  of  frequency 
fluctuations  caused  by  a  contaminant  gas  of  different  species  of  molecules  in 
a  simple  adsorption-desorption  process.  The  molecules  are  assumed  to  adsorb 
and  desorb  nondissociatively  on  the  surface  electrodes.  The  adsorbed  molecules 
are  also  assumed  not  to  exert  lateral  effects  on  neighboring  adsorbed  molecules, 
nor  do  they  diffuse  along  the  surface.  The  contaminant  gas  is  confined  within  a 
small  volume,  such  as  a  crystal  holder. 

II.  Dynamic  Equilibrium  Between  the 
Resonator  Surface  and  the  Contaminant  Gas. 

Due  to  energy  trapping,  the  mass  loading  effect  of  a  fluctuating  monolayer  of 
molecules  is  confined  almost  exclusively  to  an  area  within  the  patch  electrodes. 
The  rate  of  arrival  of  molecules  at  the  surface  can  be  evaluated  readily  from 
simple  kinetic  theory  of  gases  involving  the  pressure  and  kinetic  energy  of  the 
molecules.  If  a  site  spacing  of  0.5  nm  ;s  assumed,  the  rates  of  arrival  of  different 
molecules  at  a  site  are 

r,-  =  87.8  x  10 ^Pi/y/MiT  molecules/site/s,  (1) 

where  Pj  is  the  partial  pressure  in  torrs,  M;  is  the  molecular  weight  of  molecules 
of  type  i,  and  T  is  the  temperature  in  kelvins.  Not  all  the  arriving  molecules 
which  impinge  on  the  surface  stick  to  it.  The  rate  of  contamination  of  molecules 
of  type  i  at  a  site  is 


A,-  =  $ii'i  moiecules/site/s, 


(2) 
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where  Si  is  the  probability  of  adsorption  of  arriving  molecules  of  type  i  at 
an  uncontaminated  site.  This  probability  can  be  measured  experimentally.  A 
molecule  arriving  at  a  contaminated  site  is  assumed  to  rebound  from  the  surface. 

An  adsorbed  molecule  will  reside  on  a  site  for  a  mean  time  of  stay,  and 
subsequently  desorb  from  the  surface.  The  rate  of  desorption  per  site  is 

fii  =  Ae~’7iiT  molecules/site/s  (3) 

where  A  is  assumed  to  have  a  value  of  IxlO13,  Ei  is  the  desorption  energy  of 
molecule  of  type  i  in  kcal/mole,  and  R  is  the  gas  constant.  The  value  of  Ej  can 
be  measured  empirically. 

Our  derivation  of  the  steady-state  probabilities  of  a  site  is  essentially  that  given 
by  Langmuir,  and  discussed  in  reference  [3].  A  monolayer  or  less  of  adsorbed 
molecules  is  assumed.  When  the  system  of  contaminant  gas,  and  resonator  surface 
settles  down  to  a  steady-state  dynamic  equilibrium,  the  states  of  a  site,  which  are 
either  uncontaminated  (state  0),  or  contaminated  with  a  molecule  of  type  i  (state 
i),  can  be  described  by  steady-state  probabilities.  The  rate  at  which  a  site  enters 
state  i  is  equal  to  the  rate  at  which  it  leaves,  that  is, 

poA,-  =  pm,  (4) 

where  the  term  on  the  right-hand-side  is  the  leaving  rate.  The  coefficient  po  is  the 
probability  of  a  site  being  unoccupied,  while  pj  is  the  probability  of  a  site  being 
occupied  by  a  molecule  of  type  i.  There  are  n  types  of  contaminant  molecules  in 
the  system.  Similar!  ,  for  state  0, 

n  n 

i=i  i=i 

where  the  right-hand  term  is  the  rate  at  which  the  site  leaves  state  0.  Since  a 
site  must,  at  any  given  time,  be  in  one  of  the  states,  the  steady-state  probabilities 
must  sum  to  one,  namely, 

Y^Pi  =  L  (6) 

i=0 

Hence,  from  Eqs.(4),  (5),  and  (6),  we  obtain 

P  H 

Po  = - — 1  ™ 

n  n 

iU'+E 

;=1  ;= 1 
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where  JJ /*,-  =  fii-m,...,  nn-i -fJ-n  , 


- 2 *  (10) 

Mi 

n 

and  JJ  /xjt  =  1  if  n  =  1.  (11) 

Mi 

The  steady-state  probability  of  state  0,  po,  can  be  interpreted  as  the  fraction  of 
the  area  which  is  unconiaminated,  while  pi  is  the  fraction  of  the  area  contaminated 
with  moelcules  of  type  i. 

II.  Effects  of  Fluctuating  Pressure  on  Frequency 
Fluctuations  Caused  by  Adsorbing  and  Desorbing  Molecules. 

We  consider  the  magnitude  of  pressure  fluctuations  in  a  small  crystal  holder 
such  as  the  HC-18,  with  an  assumed  volume  and  surface  area  of  0.475  ml  and 
475  mm2,  respectively.  The  pressure  increase  due  to  one  molecule  at  300  K  is 

A  RTnA  62.359  x  300  x  6.02  x  10"23  „CI  in_17  , 

Ap  = - -  _  — - -  6.54  x  10  1  torr.  (12) 

F  V  0.000475  ' 


=  6.54  x  10-17  torr.  (12) 


If  the  adsorption  sites  are  assumed  to  have  a  linear  dimension  of  0.5  r.m, 
the  interior  surface  area  of  the  crystal  holder  would  yield  N=1.9xl015  adsorption 
sites.  At  any  given  temperature,  there  is  a  dynamic  equilibrium  in  the  number  of 
molecules  adsorbing  on,  and  desorbing  from  the  surface.  If  we  further  assume 
that  the  adsorption  and  desorption  rates  are  Ai  and  fi\,  respectively,  a  Lorentzian 
type  of  autocorrelation  function  for  pressure  fluctuations  is  obtained,  namely, 

R&p(i)  =  <r2e-(Xi+^)tf 

_2  _  AT  /  A  „\2  ^lf1!  (13) 


where  o2  =  N  (A p) 


(Ai  +  mY 
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(14) 


is  the  mean  square  pressure  fluctuation.  The  maximum  value  of  the  term 

A;/*  i 

(Ai  +  /^)2 

is  0.5,  hence,  the  maximum  root  mean  square  pressure  fluctuation  is  2.02x1 0~9 
torr.  •  If  the  ambient  pressure  is,  say,  1  millitorr,  the  root  mean  square  pressure 
fluctuation  is  still  six  orders  of  magnitude  smaller,  and  hence  would  not  be  a 
factor  in  the  resonator  frequency  fluctuations.  If  the  normal  pressure  in  the  HC- 
18  crystal  holder  at  300  K  is  as  low  as  10"6  torr,  the  root  mean  square  pressure 
fluctuation  is  still  three  orders  of  magnitude  smaller. 

The  pressure  fluctuations  would  be  larger  in  crystal  holders  smaller  than  the 
HC-18.  Consider,  for  example,  a  2  mm  diameter  crystal  holder  with  a  height  of  1 
mm,  which  we  will  refer  to  as  “HC-A”  in  our  discussions:  The  enclosed  volume 
and  surface  area  are  approximately  3.14  mm3,  and  14.1  mm2,  respectively.  Using 
(12),  the  pressure  increase  due  to  one  molecule  at  300  K  would  be  3.59xl0-13  torr. 
The  number  of  sites  is  5.64xl013.  Hence,  the  maximum  root  mean  square  pressure 
fluctuation  from  (13)  and  (14)  is  1.91x10"^  torr.  This  value  is  about  three  orders 
of  magnitude  larger  than  in  the  HC-18  crystal  holder.  The  pressure  fluctuations 
would  be  a  factor  if  the  magnitude  of  enclosed  pressure  is  of  the  order  10-6  torr. 


HI.  Pressure  in  a  Small  Crystal  Holder. 


Due  to  outgassing  and  adsorption,  the  number  of  free  flying  molecules  in  the 
confined  volume  of  a  small  crystal  holder  varies  with  temperature.  Consider  n 
types  of  molecules  with  different  molecular  weights  in  the  gas  enclosed  within 
the  crystal  holder.  The  rate  of  contamination  of  molecules  at  an  uncontaminated 
site  is  from  (1),  and  (2) 

A;  =  qiPi,  where  qi  =  87.8  x  10 ^Si/y/M^T.  (15) 


From  the  universal  gas  law,  the  partial  pressures  in  the  crystal  holder  are 

Pi  =  =  NiQ, 


V 

where  Q  = 


RT 
V  ’ 


06) 


and  Ni  is  the  number  of  moles  of  molecules  of  type  i  in  the  gas.  If  the  temperature 
is  changed  from  a  reference  temperature  To,  there  is  a  change  in  the  number  of 
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free  flying  molecules  which  is  accommodated  in  the  term  Nj.  Hence,  the  partial 
pressures  change. accordingly  to 


Pi=  (JV?  +  5  (p? -?,))«,  (17) 


where  the  superscript  0  indicates  the  state  at  reference  temperature,  and  S  is  the 
number  of  moles  of  contamination  sites.  We  subsitute  the  partial  pressures,  Pj,  of 
(15)  with  (17),  and  employing  (8),  to  yield 


A;  --  qi 


( 


iVf-i-S 


n 

Pi  ~ 


Ai  A  Hk 
kfr 


\ 


n  n  +  e  f  A;  n  n 

3= 1  V 


Q, 


(IB) 


which  are  quadratic  polynomial  functions  of  Ai.  If  there  are  n  types  of  molecules 
in  the  system,  (18)  represents  a  system  of  n  quadratic  equations  that  the  adsorption 
rates,  Aj,  must  satisfy  simultaneously.  The  system  of  equations  must  be  solved 
numerically.  However,  if  n=l,  there  is  an  exact  solution 

.  -B  ±  VB'i  -  4C 

At  —  2  > 

where  B  =  (w  +  qiQ  {S  [l  -  p?)  -  Wf})  ,  (19) 

and  C  =  fi\q\Q  (Spi  +  jV®)  . 

There  are  two  roots,  one  of  wK'ch  is  negative  and  not  meaningful.  The  positive 
root  gives  the  adsorption  rate  at  the  new  temperature  T.  Using  the  new  values  of 
Aj,  and  (15),  one  can  obtain  the  new  partial  pressures. 

We  perform  some  calculations  for  the  “HC-A”  and  HC-18  crystal  holder  of 
the  previous  section.  The  volume  of  the  crystal  holder  changes  with  temperature, 
namely, 

V  =  V0(l  +  a[T-T0)?  (20) 

where  a  is  the  coefficient  of  thermal  expansion  of  the  steel  in  the  crystal  holder, 
and  has  a  value  of  IS.SxlO^/C.  The  values  of  S  for  the  two  holders  are  9.35xl0-11, 
and  3.15xl0-9  moles,  respectively.  Some  results  of  the  pressure  of  “HC-A”  with 
one  type  of  molecules  are  shown  in  Fig.l,  where  the  molecular  weight,  sticking 
probability,  and  reference  pressure  at  25  C  are  assigned  values  of  28,  0.1,  and 
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lxlO-5  torr,  respectively.  Five  curves  are  shown,  each  exhibiting  the  pressure 
as  a  function  of  temperature  for  different  values  of  desorption  energy,  Ei.  We 
observe  that  when  the  desorption  energy  is  small,  say  Ei<9  kcal/mol,  the  system 
behaves  as  if  there  is  no  adsorption/desorption  process,  and  the  number  of  free- 
flying  molecules  remain  r  "  tively  constant  with  temperature.  Pressure  increase 
of  orders  of  magnitude,  due  to  outgassing,  is  evident  only  for  a  gas  with  larger 
desorption  energies,  say  Ei>18  kcal/mol.  Our  calculations  also  indicate  that  when 
the  initial  pressure  is  greater  than  1  torr,  the  pressure  becomes  almost  insensitive 
to  the  desorption  energy  of  contaminant  molecules. 

Figure  2  shows  the  results  for  the  HC-18  crystal  holder  with  a  gas  system 
consisting  of  five  types  of  molecules,  each  having  a  molecular  weight  of,  respec¬ 
tively,  32,  28,  16,  44,  and  18  (O2,  N2,  CH4,  CO2,  and  H2O).  Their  desorption 
energies  are  respectively  assigned  values  of  9,  12,  15,  18,  and  21  kcal/mol.  Their 
sticking  probabilities,  and  partial  pressures  at  25  C  are  all  given  a  value  of  0.1,  and 
2x10“^  torr.  The  total  pressure  at  25  C  is,  hence,  lxlO-5  torr.  We  observe  from 
Fig.2  that  the  pressure  increase  at  temperatures  above  25  C  are  principally  due  to 
outgassing  of  molecules  with  higher  desorption  energies.  The  relatively  constant 
pressure  at  temperatures  below  25  C  are  due  mostly  to  gasses  with  low  desorption 
energies.  The  composition  of  the  gas  changes  with  temperature:  At  temperatures 
above  25  C,  the  gas  has  a  higher  percentage  of  molecules  with  higher  desorption 
energies.  The  converse  is  true  at  temperatures  below  25  C.  The  results  discussed 
in  this  section  are  empirically  well  known. 

IV.  Frequency  Fluctuations  Caused  by  Adsorbing 
and  Desorbing  Molecules  in  a  Crystal  Holder. 

(a)  Continuous-Time  Markoff  Chain  for  the  Frequency  Effect  of 
a  Site,  and  Kolmogorov’s  Backward  Equations. 

The  magnitude  of  frequency  change  for  a  monolayer  of  mass  loading,  of  which 
the  mass  per  unit  area  is  much  smaller  than  the  mass  per  unit  area  of  resonator,  is 

AF;  =  fo—,  (21) 

m 

where  fo  is  the  resonator  frequency,  mj  is  the  mass  per  unit  area  of  a  monolayer 
of  molecules  of  type  i,  and  m  is  the  mass  per  unit  area  of  resonator.  Plane  waves 
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along  the  thickness  axis  were  assumed, 
the  magnitude  of.  frequency  change  per 

a{  = 

The  magnitude  of  thickness  shear  vibrations  is  maximum  at  the  center  of  resonator 
with  the  vibrations  decaying  to  zero  near  the  edges  of  electrode.  This  implies  that 
the  mass  loading  effect  of  the  monolayer  is  exerted  predominantly  by  contaminant 
molecules  located  near  the  center  of  electrode.  Hence  the  effective  Nr  value  is 
smaller  than  that  given  in  (22)  which  leads  to  larger  values  of  a*.  The  plane 
wave  assumption  is  conservative,  and  is  employed  for  the  sake  of  simplicity  in 
our  derivations. 

The  frequency  effect  fk(t)  of  a  contaminant  site  k,  in  the  process  of  adsorbing 
or  desorbing  a  molecule,  is  a  continuous-time  Markoff  chain  with  a  time  record 
given  by 


If  there  Nr  molecules  in  the  monolayer, 
site  due  to  molecule  of  type  i  is 


AF{ 
Nr  ' 


(22) 


fk  ( t )  =  ao  =  0  if  site  k  is  unoccupied,  and 
fk  ( t )  =  a;  if  site  k  is  occupied  by  a  molecule  of  type  i,  (23) 


Figure  3  shows  a  typical  time  record  of  fk(t).  The  length  of  time  a  site  remains  at  a 
certain  state  is  an  exponential  random  variable  with  a  rate  that  is  state  dependent. 
For  example,  if  site  k  is  occupied  by  a  molecule  of  type  i,  the  length  of  time  the 
site  remains  at  state  i  is  an  exponential  random  variable  with  a  rate  equal  to  the 
molecule  desorption  rate. 

We  define  a  conditional  probability  function  Py(t),  which  is  also  shown  in 
Fig.3,  as  the  probability  that  fk(t)  =  aj  given  that  fk(0)  =  a;,  that  is, 

Pij  (t)  =  P  {fk  (<)  =  aj  |  fk  (0)  =  a,}  i,j=0,l,...,n  (24) 


Since  the  molecules  adsorb  on  the  sites  independently,  the  rate  of  adsorption 
at  site  k,  uq,  is  the  sum  of  the  rates  of  adsorption  of  molecule  of  type  i,  that  is, 

n 

z/0  =  X>.  (25) 

l=i 
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The  rate  of  desorption  from  the  site,  u\,  is  dependent  on  the  type  of  molecule 
adsorbed,  namely, 


Ut  —  n jj  i-l,2,...,n. 


(26) 


The  state  of  site  k  goes  through  sequential  cycles  of  adsorption,  and  desorption 
at  rates  given  by  (25),  and  (26).  The  probability  that  an  adsorbing  molecule  is  of 
type  i  when  the  site  goes  through  a  transition  from  being  empty  to  an  occupied 
state  i  is 


VTOt  = 


At 


n  > 


j=l 


(27) 


where  7roi  represents  the  transition  probability  from  state  0  to  state  i.  The  transition 
probability  from  state  i  to  state  0  is  one 


TTjo  =  1. 


(28) 


Equations  (27),  and  (28)  cover  all  the  possible  transitions;  hence,  all  other 
transition  probabilities  are  zero 

t Uj  =  0  for  i  ^  j  ij  =  1, 2, n.  (29) 


The  conditional  probability  function  Pjj(t)  of  the  continuous-time  Markoff 
chain  can  be  shown  to  satisfy  a  system  of  first  order,  homogeneous,  differential 
equations  known  as  the  Kolmogorov’s  backward  equations[4] 


d(yPx^t^j  =  v{  KikPkj  (<)  -  ViPij  (t) ,  i,j  =  0, 1, n, 
kfr 

which,  employing  (25)  to  (29),  can  be  written  in  the  matrix  form 


where  P  = 


dP(t) 

■  =  AP 

dt 

Poo  (<) 

Poi  (t)  ■■■ 

Pon  (0 

^10  (t) 

Pn(t)  ... 

Pm  (t) 

_Pn 0  {*) 

•  •  •  #  •  • 

a  nl  ... 

Pnn  CO . 

(30) 


(31) 


(32) 
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and  A  = 


r  n 

-£At- 

Ai 

A2 

A„  ’ 

i=i 

0  ... 

“Ml 

0 

A*  2 

0 

-M2 

0 

.  Hn 

0 

0  ... 

— M«  - 

(33) 


Since  there  is  no  transition  when  the  elapsed  time  t  is  zero,  the  initial 
conditions  P(0)  for  (31)  form  an  identity  matrix.  The  Kolmogorov’s  backward 
equations  can  be  solved[5]  to  yield 

P  (t)  =  VEV-1,  (34) 


where  V  is  a  matrix  whose  column  vectors  are  the  eigenvectors  of  A  in  (33),  and 


E(*)  = 


0 

e-/M 


0  0 


0  ' 
0 

e~Pnt 


(35) 


where  terms  fax  (i=0,l,2,...,n)  are  the  eigenvalues  of  A.  Since  A  is  singular,  a 
zero  eigenvalue  always  exists,  and  we  can  define  fa  equal  to  zero.  The  matrix 
of  conditional  probabilities  P(t)  are  funcdons  of  elapsed  time  t  arising  from  the 
exponential  time  functions  in  E(t). 


(b)  Autocorrelation  Function,  and  Spectral  Density  Function 
of  Resonator  Frequency  Fluctuations. 

The  autocorrelation  function  of  the  continuous-time  Markoff  chain  f^(t)  in 
(23)  is  defined  as  the  expected  value  of  the  product  of  fk(t),  and  fk(0),  that  is, 


%(<)  =  «*  [A  »/*(<>)]■ 

Using  the  definition  of  an  expected  value[4],  we  can  write  (36)  as 


n  n 


Rfk  00  =  53  S  a'aiP  (0  =  °i»  /*  (°)  =  G»'}  • 
*'=o  j=o 


(36) 


(37) 


where  the  term  P{fk(t)=aj,  fk(0)=a; }  is  the  joint  probability  function  of  fk(t) 
and  fk(Q).  The  joint  probability  function  can  be  calculated  by  employing  the 


10 


conditional  probabilities  in  (34),  and  steady-state  probabilities  in  (7)  and  (8), 
namely, 

p{fk(t)  =  »,-,/*(  0)  =  ai}  =  (38) 

Hence,  (37)  can  be  written  as 

n  n 

%  (0  =  S  5Z  aiajPijPi-  (39) 

i=0  j=0 

The  resonator  frequency  change  is  the  algebraic  sum  of  the  frequency  effects 
of  all  contaminant  sites  on  the  two  surfaces  of  quartz  plate,  which  can  be  expressed 
as 

2  Nr 

j-m-Eaw.  (40) 

*=1 

The  autocorrelation  function  of  the  resonator  frequency  change  is,  by  definition, 
the  expected  value  of  the  product  of  F(t)  and  F(0),  that  is, 

RF(t)  =  E[F(t)F(0)}.  (41) 

When  F(t)  and  F(0)  in  (41)  are  substituted  by  the  relation  in  (40),  we  obtain 

'2  Nr  2Nr 

RF(t)  =  E  ;>>  (()£>(«)  •  (42) 

-fc=l  /=! 

Since  the  effects  of  contaminant  sites  are  mutually  independent,  the  expected 
values  of  cross  product  terms  are  zero;  and  (42)  is  simplified  to 

2Nr 

RF(t)  =  Y,Et[fk(t)fk(0)).  (43) 

*=1 

The  coefficients  on  the  right-hand  side  of  (43)  are  autocorrelation  functions  of 
frequency  effects  due  to  site  k;  therefore,  using  the  expression  in  (36), 

2  Nr 

RF(t)  =  Y,Rh(t)-  (44) 

*=1 
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We  note  that  the  generic  site  k  is  indistinguishable  from  other  sites  on  the  surfaces, 
hence,  we  can  further  simplify  (44)  to 


Rp  ( t )  =  2 NrRfk  ( t ) . 


(45) 


The  relation  for  the  autocorrelation  function  of  site  k  is  given  in  (39),  hence, 

n  n 

Rf  (t)  =  2 Nr  Y2  aiajPij  (t)  Pi-  (46) 

i=0  j= 0 

We  substitute  the  term  Py(t)  by  (34),  and  rewrite  (46)  as 

Rr(t)  =  (47) 

;=o 

n  n 

where  cr?  =  Y2  ^NrakaivkiwilPk  (48) 

k= 0  /=0 

The  coefficients  v^,  and  w^  are  the  matrix  elements  of  V,  and  V"1,  respectively. 
We  observe  that  the  autocorrelation  function  in  (47)  is  made  up  of  a  sum  of 
(n+1)  Lorentzian  autocorrelation  functions.  Each  site  has  (n+1)  states.  We 
can  now  see  that  the  reciprocals  of  eigenvalues  /?;  are  the  correlation  times 
of  the  autocorrelation  function.  The  terms  crj2  are  the  mean  square  frequency 
fluctuations,  except  for  era2  which  is  associated  with  the  zero  eigenvalue,  /?o-  The 
term  cr02  represents  a  constant  dc  term  in  the  noise  process,  and  is  2Nr  times  the 
square  of  mean  frequency  fluctuations  of  one  site.  When  it  is  removed  from  Rp(t) 
in  (47),  we  obtain  the  autocorrelation  function  of  frequency  fluctuations  about  a 
mean  frequency, 


fCO  =  ^  (49) 

;=i 

We  substitute  (22)  into  (48),  noting  that  ao=0,  to  obtain  the  mean  square  frequency 
fluctuations  in  terms  of  AFj,  namely, 


*.?  =  EE» 

Jk=l  /=! 


AFkAFt 
Nr 


■VkiWilPk 


(50) 
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Most  noise  measurements  are  made  in  terms  of  the  spectral  density  of 
frequency  or  phase  fluctuations.  In  the  present  work,  we  take  the  Fourier  transform 
of  the  autocorrelation  function  in  (49),  and  use  Wiener-Khintchine  relation  to 
obtain  the  one-sided  spectral  density  of  frequency  fluctuations  SAF(f): 


q  /  f  \  _  Y'  4ai  & 

AF{f)~hw+47r2f2Y 


(51) 


where  f  is  the  Fourier  frequency.  (49),  and  (51)  reveal  all  the  essential  infor¬ 
mation  about  the  resonator  frequency  fluctuations  induced  by  the  adsorption  and 
desorption  of  n  species  of  molecules.  They  are  similar  to  (22),  and  (23)  of  ref¬ 
erence  [2],  which  were  derived  for  one  species  of  molecules.  The  magnitude  of 
spectral  density  is  proportional  to  the  square  of  AFi,  which  is  the  mass  loading 
effect  of  a  monolayer  of  molecules  of  type  i;  and  inversely  proportional  to  Nr. 
The  magnitude  of  Nr  depends  on  the  size  of  electrode  patch,  and  contaminant 
site  spacing.  The  spectral  density  function  given  by  (51)  is  a  superposition  of  n 
number  of  Lorentzian  spectra  with  a  distribution  of  comer  frequencies  /%  given 
by  the  eigenvalues  of  (33).  Our  model  has  characteristics  similar  to  that  given 
by  McWhorter[6].  Over  a  certain  range  of  Fourier  frequencies,  1/f  noise  can 
be  generated  if  /%  is  uniformly  distributed,  having  at  least  one  comer  frequency 
per  decade  of  Fourier  frequency [7],  and  an  associated  variance,  of,  which  is  of 
significant  magnitude. 

The  spectral  density  of  frequency  fluctuations  at  1  Hz  is  useful  for  comparing 
data  of  noise  magnitudes  resulting  from  different  initial  pressures,  temperatures, 
and  types  of  molecules,  and  is  obtained  from  (51)  by  setting  f  equal  to  one, 


'S'af(I)  = 


4  of  A 

U  (/?  +  ) ' 


(52) 


Y.  Noise  Calculations  for  a  525  MHz,  UHF  AT-cut  Quartz  Resonator. 

The  fundamental  frequency  limitations  imposed  by  conventional  lapping  and 
polishing  methods  can  be  overcome  by  the  application  of  chemical  polishing 
techniques[8].  The  feasibility  of  producing  above-500  MHz  fundamental  mode 
resonators  has  been  demonstrated^].  We  consider  a  525  MHz  AT-cut  resonator 
with  an  electrode  area  of  0.003mm2,  and  thickness  of  3.2  fim.  The  resonator  is 
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confined  within  the  “HC-A”  crystal  holder,  and  pressure  fluctuations  are  neglected. 
If  we  assume  the  contaminant  sites  on  the  plate  electrodes  have  a  linear  dimension 
of  0.5  nm,  then  the  value  of  Nr  is  12xl09  sites.  The  frequency  change  for  a  0.5 
nm  monolayer  of  molecules  with  a  density  equal  to  that  of  quartz  is,  from  (21), 

m! 

AF  =  /o—  =  82  KHz,  (53) 

m 

where  fo  is  525  MHz.  Since  the  values  of  molecular  weights  to  be  considered 
are  all  of  the  same  order  of  magnitude,  we  simplify  the  calculations  by  setting 
AFi  (i=l,2,3,4,5)  to  that  given  in  (53). 

Figure  4  shows  the  spectral  density  of  frequency  fluctuations  at  1  Hz  as  a 
function  of  temperature  for  a  contaminant  gas  with  one  type  of  molecules.  Four 
curves  are  exhibited  for  desorption  energies  ranging  from  9  to  18  kcal/mole. 
The  pressure  in  crystal  holder  as  a  function  of  temperature  is  obtained  from  the 
pressures  curves  of  Fig.l.  The  initial  pressure  at  25  C  is  lxlO-5  torn  We  observe 
that  in  general  the  noise  levels  are  higher  for  molecules  with  higher  desorption 
energies.  For  temperatures  above  30  C,  a  gas  having  a  desorption  energy  of  18 
kcal/mol.  exhibits  spectral  densities  of  frequency  fluctuations  at  1  Hz  which  are 
more  than  ten  orders  of  magnitude  greater  one  with  a  desorption  energy  of  9 
kcal/mol. 

A  series  of  simulations  is  performed  to  study  the  characteristics  of  spectral 
density  of  frequency  fluctuations  given  by  (51).  Spencer  and  Smith[10]  performed 
a  spectrometric  analysis  of  residual  gases  in  metal  enclosures,  and  found  predom¬ 
inant  gases  such  as  nitrogen,  carbon  dioxide,  and  oxygen,  and  traces  of  methane, 
and  water. 

Values  of  the  sticking  probabilities,  and  desorption  energies  are  needed  in 
order  to  calculate  the  rates  of  adsorption,  and  desorption.  Accurate  values  of 
desorption  energies  of  contaminant  gases  in  quartz  resonators  are  generally  not 
available.  Water  has  desorption  energies  between  22,  and  24  kcal/mole  on  metal 
surfaces[ll].  The  sticking  probability  is  sensidvely  dependent  on  the  gas-metal 
system,  the  values  of  which  span  the  range  10-7  <  Si  <  1.  For  hydrogen, 
oxygen,  carbon  monoxide,  and  nitrogen  molecules  on  metals,  and  at  substrate 
temperatures  where  chemisorption  does  occur,  the  sticking  probabilities  [12]  are 
generally  between  0.1,  and  1.  We  assume  for  our  simulations  a  sticking  probability 
of  0.1,  and  desorption  energies  of  14.0,  15.5,  17.0,  18.5,  and  20.0  kcal/mole  for 
oxygen,  nitrogen,  methane,  carbon  dioxide,  and  water  molecules,  respectively. 
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There  are  no  compelling  reasons  for  choosing  these  values  other  than  for  the 
purpose  of  estimating  the  magnitude,  and  characteristics  of  the  noise  process.  The 
following  table  shows  the  total  pressures,  and  partial  pressures  for  six  simulations 
numbered  SI  to  S6.  The  percent  total  pressure  for  simulations  SI,  S2,  and  S3 
are  chosen  from  Spencer  and  Smith’s[8]  data  on  one  crystal  unit.  Simulations  SI 
to  S3  will  show  the  effects  of  total  pressure  of  a  contaminant  gas  mixture  with  a 
fixed  composition.  Simulations  S4  to  S6  attempt  to  show  the  effects  of  varying 
composition  of  the  contaminant  gas  mixture. 


Simulation 

number 

Total 

pressure 

at  25  C 
(torr) 

%  total 
pressure 

Oxygen 

Nitrogen 

Methane 

Carbon 

dioxide 

Water 

SI 

IE-5 

5.3 

77.7 

0.7 

15.6 

0.7 

S2 

IE-3 

5.3 

77.7 

0.7 

15.6 

0.7 

■i 

IE-1 

5.3 

77.7 

0.7 

15.6 

0.7 

S4 

IE-3 

20.0 

20.0 

20.0 

20.0 

20.0 

S5 

AE-3 

95.0 

1.25 

1.25 

1.25 

1.25 

S6 

IE- 3 

1.25 

1.25 

1.25 

1.25 

95.0 

Figure  5  exhibits  the  spectral  density  of  frequency  fluctuations  curves  for 
simulations  SI,  S2,  and  S3.  Y/e  observe  that  the  frequency  fluctuations  do  not 
follow  a  simple  power  law  noise  process:  the  slope  of  spectral  density  curves 
range  from  l/f°  to  1/f2.  The  difference  in  total  pressure  between  simulations 
SI,  and  S3  is  four  orders  of  magnitude,  while  their  spectral  densities  generally 
differ  by  about  an  order  of  magnitude.  In  Fig.6,  we  observe  that  by  changing  the 
composition  of  the  gas  according  to  the  values  given  in  the  table,  which  are  all 
within  the  same  order  of  magnitude,  the  spectral  densities  can  change  as  much  as 
two  orders  of  magnitude.  Hence,  it  seems  that  the  noise  process  is  more  sensitive 
to  the  composition  of  contaminant  gas  than  to  total  pressure. 
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ABSTRACT. 

Effects  of  multilayer  contamination  on  mean  resonant 
frequency  and  frequency  fluctuations  in  thickness-shear 
resonators  arc  studied.  A  model  based  on  mass-loading 
of  contaminant  molecules  with  adsorption  and  desorption 
rates  is  developed.  Equations  relating  change  in  mean  fre¬ 
quency  and  frequency  fluctuations  to  adsorption  and  des¬ 
orption  rates  are  derived.  Since  the  adsorption  and  desorp¬ 
tion  rates  are  functions  of  pressure  and  temperature,  change 
in  mean  frequency  and  spectral  density  of  frequency  fluctu¬ 
ations  are  studied  with  respect  to  pressure  and  temperature. 
Calculations  are  performed  for  a  10  MHz  thickness-shear 
resonator.  Frequency-temperature  and  frequency-pressure 
curves  are  plotted  for  the  10  MHz  resonator.  The  curves 
do  not  follow  a  cubic  polynomial  function  and  have  a 
magnitude  in  the  range  of  10  ppm.  The  mean  square  of 
frequency  fluctuations  under  multilayer  contamination  is 
significantly  greater  than  that  under  monolayer  contami¬ 
nation.  The  spectral  density  of  frequency  fluctuations  at 
1  Hz  is  quite  constant  in  a  wide  range  of  temperatures 
(-50  to  100°C)  when  the  values  of  heat  of  adsorption  for 
the  second  and  subsequent  layers  is  close  to  that  for  the 
first  layer.  The  magnitude  of  spectral  density  of  frequency 
fluctuations  is  about  -120  dBc  (Hz2/Hz). 


I.  Introduction. 


The  t  'ct  on  resonant  frequency  due  to  a  very  thin 
mass  layer,  such  as  one  or  more  layers  of  molecules,  ad¬ 
sorbed  on  a  thickness-shear  quartz  resonator  is  measur¬ 
able,  and  is  employed  in  quartz  microbalances  for  surface 
adsorption  studies!!].  Since  the  adsorbed  molecules  are 
known  to  have  finite  lifetimes  on  the  surface,  the  num¬ 
ber  of  adsorbed  molecules  at  a  given  instant  will  fluctu¬ 
ate.  Hence,  short-term  frequency  instabilities  will  result 
if  the  resonator  is  sufficiently  sensitive  to  the  fluctuations. 
This  was  shown  to  be  theoretically  possible  for  monolayer 
contamination  in  ultra-high-frequency  thickness-shear  res- 
onators[2].  We  examine  in  this  paper  tire  effects  of  tem¬ 
perature  and  pressure  on  changes  in  mean  resonant  fre¬ 
quency  and  frequency  fluctuations  due  to  multilayer  ad¬ 


sorption  of  contaminant  molecules  on  the  surfaces  of  a  10 
MHz  thickness-shear  resonator. 

II.  Adsorption  of  Gases  and 
Vapors  on  a  Resonator  Surface. 

For  multilayer  adsorption  of  gases  and  vapors  on 
the  surfaces  of  a  thickness-shear  resonator,  we  employ  a 
BET  model  which  was  proposed  by  Brunauer,  Emmett, 
and  Teller[3].  The  reader  is  also  referred  to  a  text  by 
Adamson[4].  The  assumptions  are:  (1)  The  number  of 
molecules  adsorbing  on  a  given  surface  is  equal  the  number 
of  molecules  desorbing  during  a  steady  state  equilibrium, 
(2)  the  heat  of  adsorption  for  the  first  layer,  £/,  has 
some  special  value,  whereas  for  all  succeeding  layers,  it  is 
equal  to  the  heat  of  condensation  of  the  liquid  adsorbate, 
Ev,  and  (3)  the  adsorption  and  desorption  can  occur  only 
from  or  on  exposed  surfaces.  Figure  1  from  reference 
[4]  shows  portions  of  surfaces  S,-  covered  by  i  (i  =  0, 
1,  2,...,)  layers  of  contaminant  molecules.  The  adsorbed 
molecules  are  assumed  to  have  no  lateral  interactions  with 
their  neighbors. 

a)  Adsorption  and  Desorption  Rates 
of  Contaminant  Molecules. 

The  rate  of  arrival  of  molecules  at  a  surface  can 
be  evaluated  readily  from  simple  kinetic  theory  of  gases 
involving  the  pressure  and  kinetic  energy  of  the  molecules. 
If  a  site  spacing  of  0.5  nm  is  assumed,  the  rate  of  arrival 
of  molecules  at  a  site  is 

r  =  87.8x10 6P/\/MT  molecules/site/s,  (1) 

where  P  is  pressure  in  tons,  M  is  molecular  w.eight  of 
molecules,  and  T  is  temperature  in  kelvins.  Not  all  the 
arriving  molecules  which  impinge  on  the  surface  stick  to 
it  The  rate  of  contamination  of  molecules  at  a  site  is 

A,-  =  s;r  molecuies/site/s,  (2) 

where  Si  is  the  sticking  coefficient  which  the  probability 
of  adsorption  of  a  molecule  impinging  at  the  surface  5',. 
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An  adsorbed  molecule  will  reside  on  a  site  for  a  mean 
time  of  stay,  and  subsequently  desorb  from  the  surface. 
The  rate  of  desorption  per  site  from  surface  Si  is 

Pi  =  molecules/site/s,  (3) 

where  A  is  assumed  to  have  a  value  of  lxlO13,  Ej  is  the 
heat  of  adsorption  in  kcal/mole,  and  R  is  the  gas  constant. 
For  surfaces  S,  where  i  is  greater  than  one,  the  rate  of 
desorption  per  site  is 

/i2  =  Ac~e'!'^  molecules/site/s,  (4) 


III.  Changes  in  Mean 
Resonant  Frequency. 

From  Eq.(7),  we  observe  that  the  steady-state  proba¬ 
bilities  pi,  being  functions  of  the  adsorption  and  desorption 
rates,  are  functions  of  temperature  and  pressure  for  a  given 
absorbate-contaminant  system.  The  sensitivity  of  the  res¬ 
onator  to  mass  loading  is  confined  mainly  to  the  electrode 
patch.  A  monolayer  of  contaminant  molecules  covering 
the  electrode  patch  will  change  the  resonant  frequency  by 
the  amount 


where  Ev  is  the  heat  of  condensation  in  kcal/mole. 


b)  Dynamic  Equilibrium  Between  the 
Resonator  Surface  and  the  Contaminant  Gas. 

The  derivation  of  the  steady-state  probabilities  of  a 
site  being  uncontaminated,  or  contaminated  with  one  or 
more  molecules,  follows  essentially  that  given  by  the  BET 
model.  We  presently  assume  that  the  adsorbed  mass  layer 
will  consists  of  no  more  than  three  layers  of  molecules. 
Hence,  the  sticking  coefficient  sj  is  given  a  value  of 
zero.  Derivations  for  models  with  more  than  three  layers 
are  similar  to  what  is  given  here.  When  the  system  of 
contaminant  gas  and  resonator  surface  settles  down  to  a 
steady  state  dynamic  equilibrium,  the  slates  of  a  site,  which 
is  either  uncontaminated  (state  0),  or  contaminated  with 
i  molecules  in  a  layerwise  manner,  can  be  described  by 
steady-state  probabilities  p,-  (i  -0,1,2,  3).  The  rate  at 
which  a  site  enters  state  i  is  equal  to  the  rate  at  which  it 
leaves,  that  is, 

PoA0  =  P1P1 

PlAl=P2P2  (5) 

P2A2  =  P3P2- 

The  steady-state  probabilities  must  sum  to  one,  since  the 
site  must  at  any  given  time  be  in  one  of  the  states, 


Po  +  pi  +  P2  +  P3  —  1- 


(6) 


The  steady-state  probabilities  p,-  can  be  expressed  in  terms 
of  the  adsorption  and  desorption  rates  using  Eqs.(5)  and 
<6): 

n=D 
Ao 


Pi  = 


P  2  = 


AqAi 
P1P2  D 
A0A1A2 


(7) 


P3  pipP 

where  z?=i  +  ^  +  Mi  +  Ml^. 

Pi  P1P2  P1P2 


,  „  m' 

Ay  =  -/o —  (8) 

m 

where  m’  and  m  are  the  mass  per  unit  area  of  contami¬ 
nant  molecules  and  mass  per  unit  area  of  resonator  plate, 
respectively.  The  resonator  frequency  is  fo.  The  negative 
sign  signifies  that  the  resonant  frequency  decreases  with 
mass  loading.  Eq.(8)  is  valid  only  when  the  mass  per  unit 
area  of  contaminant  molecules  is  much  smaller  than  the 
mass  per  unit  area  of  resonator  plate.  The  resonator  is 
assumed  to  vibrate  in  a  pure  thickness-shear  mode.  The 
actual  change  in  frequency  will  depend  on  the  number  of 
molecular  layers  and  fraction  of  surface  covered,  namely, 

A/  = 

fo  (9) 

2Ay  ((pi  -  Plo)  +  2  (p2  -  P20)  +  3  (p3  -  p3o)) 

where  p,0  (i  =  1 , 2, 3)  are  the  steady-state  probabilities  at 
the  reference  temperature  and  pressure,  and  the  L;.  u  :  ; 
on  the  right  hand  side  accounts  for  containing  ■  ha. 
top  and  bottom  surface. 

a)  Changes  in  Mean  Resonant  Frequency  for 
a  10  MHz  Thickness-Shear  Resonator. 

Simple  calculations  are  performed  for  a  10  MHz 
thickness-shear  resonator  to  examine  the  characteristics 
and  magnitude  of  changes  in  mean  frequency  as  a  func¬ 
tion  of  temperature  and  pressure.  The  resonator  thickness 
is  0.165  mm  and  the  area  of  circular  electrode  patch  is  10 
mm2.  A  0.5  nm  thick  monolayer  of  molecules  covering 
the  top  electrode  and  having  a  density  equivalent  to  quartz 
will  yield 

Ay  =  30  Hz, 

Ay  „  (10) 

or  — ' L  =  3  ppm. 
fo 

The  molecular  weight  of  contaminant  molecules  is  taken 
as  equal  to  28.  For  our  calculations  of  the  adsorption  rates, 
we  take  values  of  sticking  coefficients  s;  (i  =  0,1,2)  equal 
to  0.1. 
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Figure  2  shows  the  effects  of  surface  contamina¬ 
tion  with  temperature  on  the  mean  resonator  frequency 
in  parts  per  million.  Three  frequency-temperature  curves 
are  given  corresponding  to  heats  of  adsorption  15, 18,  and 
21  kcal/mole.  The  heats  of  condensation  are  taken  to  be 
75  percent  of  their  respective  heats  of  adsorption.  The 
pressure  is  kept  constant  at  0.001  torr.  We  observe  that 
the  surface  contamination  as  a  function  of  temperature  can 
cause  the  mean  frequency  to  change  as  much  as  10  ppm. 
The  f-T  curves  do  not  follow  a  cubic  polynomial  func¬ 
tion.  The  characteristics  and  magnitudes  of  curves  will  be 
affected  by  initial  pressures,  heats  of  adsorption  and  con¬ 
densation,  and  sticking  coefficients.  The  mean  frequency 
increases  with  increasing  temperature. 

Figure  3  shows  the  effects  of  surface  contamination 
with  pressure  on  the  mean  resonator  frequency  in  parts 
per  million.  The  three  frequency-pressure  curves  have  the 
same  corresponding  heats  of  adsorption  and  condensation 
as  in  Fig.2.  The  temperature  is  kept  constant  at  25°C. 
We  observe  that  the  mean  frequency  is  also  sensitive  to 
pressure  changes,  and  is  of  a  higher  order  function  than 
a  cubic  polynomial.  The  mean  frequency  decreases  with 
increasing  pressure. 


IV.  Theoretical  Developments 
for  Frequency  Fluctuations. 


Although  there  is  a  steady-state  equilibrium  in  the 
adsorption  and  desorption  of  contaminant  molecules,  there 
are  instantaneous  fluctuations  in  the  number  of  adsorbed 
molecules.  Each  contaminant  site  on  the  electrode  patch  of 
a  thickness-shear  resonator  exerts  a  small  but  finite  mass 
loading  effect  on  the  resonant  frequency.  Hence,  the  fluc¬ 
tuations  will  be  manifested  in  short  term  frequency  insta¬ 
bilities.  We  investigate  the  magnitude  and  spectral  char¬ 
acteristics  of  short  term  frequency  fluctuations  as  caused 
by  spatial  and  temporal  fluctuations  in  the  number  of  ad¬ 
sorbed  contaminant  molecules. 


a)  Birth  and  Death  Processes. 


Birth  rate  in  So  sites  =  Ao 
Birth  rate  in  Si  sites  =  Aj 
Birth  rate  in  52  sites  =  A2 
Death  rate  in  Si  sites  r 
Death  rate  in  S2  sites  =  r 
Death  rate  in  S3  sites  =  /i2 

The  birth  and  death  process  is  a  continuous-time 
Markov  chain  which  is  more  generally  an  exponential 
model.  Hence,  the  state  transition  rate  for  a  5;  site  is 
simply  the  sum  of  birth  and  death  rates  associated  with 
the  state  of  the  site,  that  is, 


50  site:  A0 

51  site:  Ai  +  /xj 

52  site:  A2  +  /x2 

53  site: 


(12) 


b)  Frequency  Effect  of  a  sit e/*(7). 

Since  each  contamination  site,  say  site  k,  exerts  a 
finite  effect  on  the  resonant  frequency,  the  frequency  effect 
can  be  described  by  a  continuous  function  fk(t)  which  is  a 
continuous-time  Markov  chain.  The  function  has  specific 
values  depending  on  the  state  of  site: 

fk  ( t )  =  a0  =  0  for  S0  site, 

fk  ( t )  =  ai  =  jf-  for  Si  site, 

fk  (t)  =  a2  =  for  52  site,  (1'  ; 

yvr 

fk  ( t )  =  a3  =  for  S3  site. 

where  N,  is  the  number  of  sites  on  one  electrode  patch. 
The  magnitude  of  thickness-shear  vibrations  is  maximum 
at  the  center  of  electrode  patch  with  the  vibrations  de¬ 
caying  to  zero  near  the  edges  of  electrode.  A  graph  of 
the  magnitude  of  thickness-shear  vibrations  will  exhibit 
a  Gaussian  profile.  This  implies  that  the  mass  loading  ef¬ 
fect  of  the  contaminant  molecules  is  exerted  predominantly 
near  the  center  third  of  electrode.  Hence  the  effective  Nr 
value  is  smaller  than  the  number  of  sites  in  one  electrode 
patch.  Subsequent  analysis  will  show  that  a  smaller  Nr 
value  leads  to  larger  magnitudes  of  root-mean-square  fre¬ 
quency  fluctuations. 


The  stochastic  analysis  of  fluctuations  in  the  ad¬ 
sorption  and  desorption  of  contaminant  molecules  fol¬ 
lows  the  usual  developments  for  birth  and  death  pro- 
cesses[5].  The  table  below  gives  the  birth(adsorption)  and 
death(desorption)  rates  which  follow  from  Eqs.(2)  to  (4): 


c)  Conditional  Probability  Function 
Pip)  of  fk(t). 

The  probability  of  fk(t)  at  a  subsequent  time  1  having 
a  value  of  aj,  given  an  initial  value  of  a,  at  t  equal  to  zero, 
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« 


is  defined  by  a  conditional  probability  function 
Pij(t)  =  P{h(t)  =  aj\fk(0)  =  ai} 
i>j  =  0)  1|2,3. 


(14) 


For  a  continuous-time  Markov  chain,  Eq.(14)  must  satisfy 
Kolmogorov’s  backward  equations  whicli  are  a  system  of 
first  order  differential  equations: 


(15) 


■  ><') 

=  AP((), 

where  P  (l)  = 

rPooW  Poi(t) 

P02(t) 

P03  (0l 

Pio(()  Pn  (t) 

Pn(t) 

Pl3(0 

P20{t)  Putt) 

P22W 

P23  (t) 

lP3o(t)  Pai  (() 

P32W 

P33  (<)  J 

and  A 

.  = 

© 

o 

< 

1 

1 _ 

0 

0 

Pi 

-(Ai  +  p\) 

Ai 

0 

0 

P2 

-(A2  +  /12)  ^2 

.  0 

0 

P-2 

~P2 

(16) 


(17) 


The  initial  conditions  for  Eq.(15)  is  an  identity  matrix 
since  at  zero  time  there  are  no  transitions. 


P(0)  = 


1 

0 

0 

LO 


01 

0 

0 

1J 


(18) 


d)  Autocorrelation  Function  of 
Frequency  Fluctuations  about  Mean 
Resonant  Frequency. 

The  autocorrelation  function  of  the  continuous-time 
Maikov  chain  fk(t)  in  £q.(13)  is  defined  as  the  expected 
value  of  the  product  of  fk(t),  and  fk(0),  that  is, 

%(*)  =  £  [/*(*)/*  (0)1,  (22) 

which  by  definition  can  be  written  as 

R-h  (t)  = 

3  3  (231 

Y  Y a‘aiP  Uk  w  =  ai ' h  (°)  =  a<) ' 

i=0j=0 

where  the  term  P {fk  (t)  =  aj,fk  (0)  =  a,}  is  the  joint 
probability  distribution  function  of  fk(t),  and  fk(0).  The 
joint  probability  distribution  function  can  be  calculated 
by  employing  the  conditional  probabilities  in  Eq.(19),  and 
steady-state  probabilities  in  Eq.(7).  Hence, 


P  ih  (0  =  a; ,  ft  (0)  =  a,}  =  P\j  ( t )  Pi-  (24) 


Using  Eqs.(24)  and  (21),  Eq.(23)  can  be  written  as 

%  (0  =  Y  rte~Unt  (25> 

m= 0 


The  Kolmogorov’s  backward  equations  can  be  3  3 

solved[6)  to  yield  where  <rm=)  ,  /  yaiajVimVmjPi.  (26) 


where 

P(f)  =  VEV"1 

-e-ai°j  0  0 

0  p-vi*  0 

0  • 

0 

(19) 

t=0 j= 0 

Since  wo  is  the  zero  eigenvalue,  is  a  constant 

dc  noise  which  is  the  square  of  mean  value  of  fk(t).  In 
our  study  of  short  term  instabilities  in  frequency,  we  are 
interested  in  frequency  fluctuations  about  a  mean  resonant 

E(f)  = 

V  c  u 

0  0  e~W3t 

.  0  0  0 

V/ 

0 

e-a’3<. 

(20) 

frequency.  Hence,  we  define  a  continuous-time  Markov 
chain  of  frequency  fluctuations 

is  a  diagonal  matrix  of  exponential  time  functions  with 
constants  wm  (m-0, 1,2,3)  which  are  the  absolute  values 
of  eigenvalues  of  A  in  Eq.(17).  Since  A  is  singular,  a 
zero  eigenvalue  exists,  and  we  can  define  wo  to  be  the 
zero  eigenvalue.  All  the  nonzero  eigenvalues  are  negative. 
The  column  vectors  vm  of  V  and  eigenvalues  -w„  form 
the  eigenpairs  {-w„,  of  A.  Eq.(19)  and  (20)  show 
that  the  matrix  of  conditional  probabilities  are  exponential 
functions  of  elapsed  time  t,  namely, 

Pii  (<)  =  E  (21) 

m= 0 

where  v,m  and  t»~j  are  elements  of  matrices  V  and  V-1, 
respectively. 


A  fk(l)  =  fk(t)-c0  (27) 

The  autocorrelation  function  of  frequency  fluctuations 
about  a  mean  frequency  for  site  k  is 

R*/k(t)  =  E[{Afh(t)}{Afk(0)}} 

=  E[fk(t)fk(0)]-<rl  .. 

3  M 

=  Y 

m= 1 

Note  that  since  ao  is  equal  to  zero,  Eq.(26)  can  also  be 
written  as 

=  YY  a'ai  V*mVmjPi  ■  (29) 

»=1  j=l 
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We  observe  that  the  sum  of  <r over  m  =  1,2,3  is  the 
mean  square  of  frequency  fluctuations  due  to  site  k. 

There  are  2 Nr  sites  in  the  resonator,  which  accounts 
for  contamination  sites  on  the  top  and  bottom  electrodes. 
Hence,  the  resonator  frequency  fluctuation  about  a  mean 
frequency  at  any  tim°.  t  is 


2Nr 

.  A F(t)  =  X>/*(<), 

t=i 

whose  autocorrelation  function  is 
RAF(t)  =  E[AF(t)AF(0)], 

‘/2Nr  \  /2  Nr  N 

=  e  EA/<(°) 

.  \jt=i  /  Vi=i  t 

2  Nr  2  Nr 

=  EE  EUAMtMAfdO)}} 

t=i i=i 


(30) 


(31) 


Since  the  effects  of  contaminant  sites  are  mutually  inde¬ 
pendent,  the  expected  values  of  cross  product  terms  are 
zero,  and  Eq.(31)  becomes 


2ATr 

R*r(i)  =  22El{  A/fc(0HA/*(0)}] 

J*  (32) 

=  X>/*(0 

t=l 

The  generic  site  k  is  indistinguishable  from  other  sites 
on  the  electrode  surfaces,  therefore,  the  autocorrelation 
function  in  Eq.(32)  is 


R*F(i)  =  2NrRth(i).  (33) 


Hence  the  root-mean-square  of  frequency  fluctuations  is 
proportional  to  the  frequency  change  of  resonator  to  a 
monolayer  of  contaminant  molecules,  and  inversely  pro¬ 
portional  to  the  square  root  of  the  number  of  contaminant 
sites  in  electrode  patch. 


e)  One-Sided  Spectral  Density  of 
Frequency  Fluctuations  about  Mean 
Resonant  Frequency. 

Wiener-Khintchine  relation  is  employed  to  obtain  the 
one-sided  spectral  density  of  frequency  fluctuations  from 
the  autocorrelation  function  of  Eq.(34);  which  yields 


Saf(/)  = 


wm  +  4rr2/2  ’ 


(37) 


where  /  is  the  Fourier  frequency.  Eq.(37)  reveals  all  the 
essential  characteristics  of  frequency  fluctuations  induced 
by  multilayer  surface  contamination.  The  spectral  den¬ 
sity  function  is  a  superposition  of  three  Lorentzian  spectra 
with  a  distribution  of  comer  frequencies  um.  The  function 
is  similar  to  that  given  by  McWhorter[7].  Over  a  certain 
range  of  Fourier  frequencies,  1/f  noise  can  be  generated  if 
wn  is  uniformly  distributed  with  at  least  one  comer  fre¬ 
quency  per  decade  of  Fourier  frequency[8],  and  the  asso¬ 
ciated  cr having  a  significant  magnitude.  The  magnitude 
of  spectral  density  is  proportional  to  the  square  of  A/,  and 
inversely  proportional  to  Nr. 


V.  Noise  Calculations  for  a  10 
MHz  Thickness-Shear  Resonator. 


Using  the  last  equation  of  Eq.(28),  we  can  write  Eq.(33) 
as 

=  (34) 

where  a2m  -  2Nra2m. 


Noise  calculations  for  the  10  MHz  thickness-shear 
resonator  previously  defined  in  section  ni(a)  is  performed. 
In  addition,  an  average  site  spacing  of  0.5  nm  is  assumed 
for  the  contaminant  sites  on  electrode  patches.  Hence,  the 
number  of  sites  in  one  electrode  patch  is 


We  observe  that  the  autocorrelation  function  in 
Eq.(34)  is  made  up  of  ‘hree  Lorentzian  autocorrelation 
functions.  The  reciprocals  of  um  are  the  correlation  times 
of  the  autocorrelation  function.  The  term 


a 


(35) 


is  the  root-mean-square  of  resonator  frequency  fluctua¬ 
tions,  and  from  Eqs.(13)  and  (29), 


3  3  3 


E  EE2,7^vmiP.-- 

m=l  i=l  j=l 


(36) 


Nr  =  4x1 013  sites.  (38) 

Calculations  for  both  monolayer  and  multilayer  contami¬ 
nations  are  performed.  The  model  for  monolayer  contam¬ 
ination  was  previously  reported  in  reference  [2],  and  is 

pr'ltidlltt  p  oru»r>iolWoH  r%f  nwmnt 

itWkUUii^  I*  MpWiUlUjW  VUJV  v*  MIV  UtVJVIU  UIVUVIl 

Figure  4  shows  the  spectral  density  curves  of  fre¬ 
quency  fluctuations  plotted  as  a  function  of  Fourier  fre¬ 
quency.  Four  curves  are  shown.  The  solid  curve  represents 
monolayer  contamination,  while  the  three  other  curves  are 
for  multilayer  contamination  using  values  of  heat  of  con¬ 
densation  Ev  which  are  0.5, 0.75  and  1.0  of  Ej,  the  heat  of 
absorption  with  a  value  of  21  kcal/mole.  We  observe  that 
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the  magnitude  of  spectral  density  for  Ev  greater  than  0.75 
Ei  is  about  50  dBc  greater  than  the  curve  with  monolayer 
contamination.  For  the  particular  pressure  (0.001  torr), 
temperature  (25°C),  and  Fourier  frequency  range  of  less 
than  100  Hz,  the  spectral  density  curves  with  Ev  less  than 
0.5  Ei  resemble  the  monolayer  curve. 

The  spectral  density  of  frequency  fluctuations  at  1 
Hz  offset  is  useful  for  comparing  data  of  noise  magnitudes 
resulting  from  different  initial  pressures,  temperatures,  and 
heats  of  absorption  and  condensation.  This  is  obtain  from 
Eq.(37)  by  setting  /  equal  to  one, 


Saf  (1)  = 


(39) 


Figure  5  shows  the  spectral  density  of  frequency  fluc¬ 
tuations  at  1  Hz  plotted  as  a  function  of  pressure.  We 
observe  that  for  wide  ranges  of  pressures  (1CTS  to  10  torr) 
the  spectral  density  with  multilayer  contamination  and  Ev 

>  0.75 Ej  has  a  magnitude  greater  than  monolayer  contami¬ 
nation.  The  magnitude  difference  is  greater  at  higher  pres¬ 
sures.  The  curve  with  Ev  =  0.5 Ei  follows  the  monolayer 
curve  up  to  a  pressure  of  about  0.002  torr,  and  diverges 
from  it  at  higher  pressures. 

Figure  6  shows  the  spectral  density  of  frequency  fluc¬ 
tuations  at  1  Hz  plotted  as  a  function  of  temperature.  Cal¬ 
culations  show  that  for  multilayer  contamination  with  Ev 

>  0 .9Ej,  the  magnitude  of  spectral  density  at  1  Hz  re¬ 
mains  remarkably  constant  in  a  wide  range  of  temperature 
(-50  to  100°C),  as  opposed  to  the  monolayer  curve.  With 
Ey  =  0J5Ej,  the  multilayer  contamination  curve  remains 
constant  in  a  shorter  range  of  temperatures  (-50  to  35°C). 
The  curve  with  Ev  =  0.5E;  follows  the  monolayer  curve 
in  temperatures  greater  than  10°C. 


VI.  Summary. 

A  model  based  on  mass-loading  of  contaminant 
molecules  having  adsorption  and  desorption  rates  is  devel¬ 
oped.  Calculations  are  made  for  a  10  MHz  thickness-shear 
resonator.  The  model  brings  out  some  interesting  aspects 
of  contamination  in  resonators  which  may  be  tested  em¬ 
pirically: 

1.  The  change  in  surface  contamination  layers  due  to  a 
change  in  temperature  affects  the  mean  resonant  fre¬ 
quency  of  thickness-shear  resonators.  This  frequency 
change  may  be  masked  by  changes  in  material  prop¬ 
erties,  such  as  elastic  stiffnesses  and  density.  Calcula¬ 
tions  for  the  10  MHz  resonator  yield  mean  frequency 
changes  of  as  much  as  10  ppm. 


2.  The  mean  resonant  frequency  is  also  affected  in  the 
same  order  of  magnitude  by  the  ambient  pressure 
which  influences  the  characteristics  of  the  surface 
contamination  layers. 

3.  Magnitude  of  spectral  density  of  frequency  fluctua¬ 
tions  for  multilayer  adsorption  and  desorption  is  larger 
than  for  monolayer  adsorption  and  desorption.  This 
is  especially  true  when  the  value  of  Ev  is  close  to  Ej. 
Calculations  for  the  10  MHz  resonator  show  that  the 
magnitude  of  mean  square  frequency  fluctuations  is 
in  the  measurable  range  of  about  -120  dBc  (Hz2/Hz). 

4.  Spectral  density  of  frequency  fluctuations  at  1  Hz  for 
multilayer  adsorption  and  desorption,  when  the  value 
of  Ey  is  close  to  £;,  is  remarkably  constant  with 
changes  in  temperature.  This  is  not  so  for  monolayer 
adsorption  and  desorption. 
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Fig.l  Adsorption  of  gases  and  vapors  on  solid  surfacc[4], 
the  BET  model(3J. 


f  Spectral  density  of  freq.  fluctuoticns 
vs  Fourier  frequency.  Pressure=JE-3  torr 
Temperature  =  25  deg  C.  El  =  21  kcot/mol 
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